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FOREWORD 


The  work  described  in  this  Annual  Technical  Report  was  performed  at 
the  United  Technologies  Research  Center  (UTRC)  under  Contract  N00014-79-C-0344 
entitled  "Study  of  Two-Phase  Nozzle  Theory  and  Parametric  Analysis",  for  the 
Office  of  Naval  Research  (ONR).  This  report  summarizes  results  obtained  for 
the  Phase  I  -  Two-Phase  Nozzle  Theory  study  program.  Dr.  Simion  C.  Kuo  is 
the  Principal  Investigator  for  this  contract,  and  Dr.  C.  W.  Deane  is  the 
major  contributor  for  this  phase  of  the  study. 

The  contract  program  was  initiated  with  ONR  on  May  15,  1979,  and  the 
ONR  Program  Manager  is  Mr.  M.  Keith  Ellingsworth,  Power  Program,  ONR,  Arlington, 
Virginia.  Valuable  guidance  and  comments  received  from  Mr.  Ellingsworth  and 
Mr.  John  R.  Satkowski,  Director  of  Power  Program  at  ONR  are  gratefully 
acknowledged. 
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Two-Phase  Theory  and  Parametric  Analysis 
Phase  I  -  Two-Phase  Nozzle  Theory 


SUMMARY 


The  object  of  the  work  presented  in  this  report  was  to  formulate  a 
working  theory  with  which  the  performance  of  two-phase  nozzles  can  be  analyzed 
parametrically  over  a  wide  range  of  flow  conditions  and  fluid  properties. 

As  a  basis  for  formulating  the  theory,  the  technical  status  of  two-phase 
nozzles  was  assessed,  the  governing  flow  parameters  were  identified,  and  the 
flow  domains  were  classified  with  the  previously  available  theory.  A  working 
theory  of  two-phase  nozzles  suitable  for  parametric  analyses  and  optimization 
was  then  formulated,  incorporating  the  significant  interphase  effects  of 
velocity  slip,  droplet  breakup,  and  droplet  heat  transfer. 

The  technical  status  of  two-phase  nozzles  was  assessed  by  reviewing 
previous  investigations,  both  analytical  and  experimental,  and  the  ranges  of 
parameters  considered  in  previous  investigations  were  summarized.  The 
limitations  of  the  previous  models  were  documented,  and  areas  where  the 
theory  could  be  improved  or  simplified  for  greater  physical  understanding  of 
two-phase  flow  were  reviewed.  Liquid/gas  combinations  for  two-phase  nozzles, 
and  suitable  combinations  for  marine  propulsion  and  other  selected  applications 
were  discussed.  Single-phase  (or  homogeneous)  theory  was  also  reviewed  from 
the  viewpoint  of  energy  conversion  in  a  nozzle  and  was  compared  with  models 
applicable  to  separate  phases  in  terms  of  their  flow  domains. 

Finally,  a  working  theory  of  two-phase  nozzles  was  formulated  based 
on  the  governing  flow  parameters  and  on  the  six  basic  equations  that  describe 
the  two-phase  flow.  From  this  theory,  a  computer  model  was  developed  which 
includes  the  significant  interphase  effects  of  velocity  slip,  droplet  breakup, 
and  heat  transfer  and  which  is  suitable  for  parametric  analyses  and  optimization. 
Preliminary  parametric  calculations  based  on  use  of  this  model  were  compared 
with  results  from  previous  two-phase  nozzle  models  and  with  available  experi¬ 
mental  data. 

This  study  program  was  conducted  by  the  Energy  Conversion  Systems 
Analysis  Group  at  UTRC  under  Contract  N00014-79-C-0344  from  the  Office  of 
Naval  Research,  Power  Program,  Arlington,  Virginia. 
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INTRODUCTION 


The  flow  of  two-phase  mixtures  through  a  nozzle  is  a  complicated 
process  which  is  receiving  increasing  attention  since  many  engineering 
applications  involve  flows  of  a  substance  that  is  a  suspension  of  liquid 
droplets  or  solid  particles  in  a  gas.  A  typical  example  is  a  two-phase 
nozzle  where  a  gas  and  a  liquid  at  a  condition  of  high  pressure  and  low 
velocity  are  mixed  at  a  nozzle  inlet  and  then  expanded  through  the  nozzle  to 
a  condition  of  low  pressure  and  high  velocity.  Two  types  of  two-phase  flow 
are  of  interest  here.  "One-component"  flow  is  the  case  where  the  gas  phase 
is  the  vapor  of  the  liquid  being  accelerated,  and  "two-component"  flow  is  the 
case  where  the  gas  phase  is  of  a  different  chemical  species  from  that  of  the 
liquid.  Among  other  applications,  two-phase  nozzles  can  be  integrated  with 
two-phase  turbines  for  prime-mover  applications  and  for  geothermal  power 
generation . 

The  major  difference  between  single-phase  (gas)  nozzle  flow  and  two-phase 
nozzle  flow  is  the  interaction  between  the  phases:  the  transfer  of  momentum 
and  heat  between  the  two  phases.  In  two-phase  nozzle  flow,  as  the  gas 
expands,  drag  forces  transfer  momentum  from  the  gas  phase  to  accelerate  the 
liquid  droplets,  while  heat  is  being  transferred  between  the  liquid  and  the 
gas.  However,  the  temperature  of  the  gas  phase  decreases  as  the  gas  expands, 
therefore  heat  transfer  from  the  liquid  phase  will  partly  offset  this 
decrease,  so  the  temperature  decrease  of  the  gas  in  a  two-phase  flow  is 
less  than  that  for  the  pure  gas  expanded  through  the  same  pressure.  Hence, 
interphase  effects  must  be  considered  in  analyzing  two-phase  nozzle  flows. 

Previous  analytical  models  of  two-phase  nozzle  flow  generally  display 
one  of  two  drawbacks.  Either:  only  some  of  the  significant  interphase 
effects  of  velocity  slip,  droplet  breakup,  and  droplet  heat  transfer  were 
considered  (e.g.,  Refs.  1  to  5);  or  these  three  effects  were  considered,  but 
the  model  (e.g.,  Ref.  6)  was  generally  formulated  with  a  view  towards  detailed 
design  and  not  with  the  intention  of  conducting  parametric  analyses  and 
optimization  where  many  cases  must  be  quickly  examined  yet  with  reasonable 
accuracy.  Therefore,  to  avoid  these  drawbacks,  a  working  theory  should  be 
formulated  which  is  suitable  for  parametric  analyses  and  optimization  yet 
includes  the  significant  interphase  effects  of  velocity  slip,  droplet  breakup, 
and  droplet  heat  transfer  as  well  as  provides  details  of  the  local  flow 
conditions  throughout  the  nozzle. 

The  object  of  the  Phase  I  -  Two-Phase  nozzle  theory  study  program  was 
to  formulate  a  working  theory  of  two-phase  nozzle  flows  for  use  in  analyzing 
nozzle  performance  as  a  function  of  the  governing  flow  parameters;  the  results 
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of  this  study  program  are  presented  in  this  report.  The  technical  status 
of  two-phase  nozzles  was  assessed  and  summarized  in  Task  1.1,  and  the  governing 
flow  parameters  and  their  ranges  of  interest  were  identified  in  Task  1.2. 

The  flow  domains  of  two-phase  nozzles  were  classified  and  correlated  in  Task 
1.3  to  identify  the  physical  characteristics  of  these  domains.  Based  on  this 
work,  a  working  theory  of  two-phase  nozzle  flow  then  was  formulated  that 
includes  the  significant  interphase  effects  and  that  is  suitable  for  parametric 
analyses  and  optimization  while  still  providing  detailed  local  flow  conditions 
in  the  nozzle. 
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RESULTS  AND  CONCLUSIONS 


1.  A  working  theory  of  two-component  two-phase  (gas-liquid)  nozzles  was  developed 
which  includes  the  interphase  transport  of  momentum  and  thermal  energy  but  with¬ 
out  mass  transfer  and  frictional  losses.  This  working  theory  provides  reason¬ 
ably  rapid  yet  accurate  computerized  numerical  solutions  suitable  for  parametric 
and  optimization  analyses  of  these  nozzles  with  gas-liquid  flows. 

2.  Previous  two-phase  nozzle  theories  generally  display  one  of  two  drawbacks: 
either  they  considered  only  some  of  the  interphase  transport  effects  for  sim¬ 
plicity  of  solution,  or  they  used  excessively  tedious  computational  procedures 
to  calculate  design  details  of  individual  two-phase  nozzles. 

3.  The  departure  of  two-phase  nozzle  behavior  from  that  of  a  single-phase  (gas) 
nozzle  can  be  attributed  to  the  interphase  coupling  of  momentum  and  thermal 
energy  transport.  As  a  result,  liquid  droplets  are  accelerated  by  the  expanding 
gas  while  heat  transfer  from  the  droplets  retards  expansion  cooling  of  the  gas 
which  could  expand  nearly  isothermally  at  a  high  liquid-to-gas  loading  ratio 
and  with  small  droplets. 

4.  The  principal  interphase  coupling  effects  which  must  be  considered  in  formulat¬ 
ing  the  two-phase  nozzle  theory  are:  1)  droplet  drag  due  to  velocity  slip; 

2)  droplet  breakup  as  dictated  by  the  aerodynamic  and  surface  tension  forces; 

3)  heat  transfer  between  droplets  and  the  surrounding  gas;  and  4)  mass  transfer. 
The  flow  parameters  which  were  identified  as  being  significant  are:  velocity 
and  temperature  of  the  liquid  and  gas;  pressure  of  the  gas;  the  liquid-to-gas 
loading  ratio;  the  droplet  diameter;  and  the  thermophysical  properties  of  the 
two  phases. 

5.  Use  of  a  hot  liquid  at  a  temperature  substantially  higher  than  that  of  the  gas 
will  increase  the  enthalpy  of  the  two-phase  mixture  as  well  as  heat  transfer 
from  the  droplets  to  the  surrounding  gas.  Additionally,  an  external  heater 
required  to  heat  the  liquid  would  be  relatively  small.  Therefore,  the  potential 
application  of  such  a  nozzle  as  an  integral  component  of  a  two-phase  turbine 
would  be  enhanced. 

6.  Among  the  various  liquid-gas  combinations  reviewed,  a  high-temperature  synthetic 
heat  transfer  oil  (such  as  one  of  those  in  the  Dowtherm  family)  suspended  in 
steam  is  an  attractive  potential  combination  for  two-phase  turbine  applications. 
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Two-Phase  Nozzle  Theory  and  Parametric  Analysis 
Phase  I  ~  Two-Phase  Nozzle  Theory 


The  overall  objective  of  this  study  program  is  to  advance  the  basic 
understanding  of  two-phase  nozzles.  To  meet  this  overall  objective,  the 
specific  objective  of  Phase  I  is  to  develop  a  working  theory  of  two-phase 
nozzles,  based  on  the  principal  governing  flow  parameters  and  on  previous  work 
on  two-phase  nozzles  in  the  literature.  Phase  I  consists  of  four  tasks:  1.1) 
Technical  Status  of  Two-Phase  Nozzles;  1.2)  Identification  of  Governing  Flow 
Parameters;  1.3)  Classification  of  Flow  Domains  with  Available  Theory;  and  1.4) 
Formulation  of  Working  Theory  for  Two-Phase  Nozzles.  The  results  obtained 
in  each  of  these  four  tasks  are  presented  in  this  report. 


1.1  Technical  Status  of  Two-Phase  Nozzles 

A  nozzle  is  a  flow  channel  in  which  the  velocity  of  the  fluid  is  increased 
and  the  pressure  is  reduced,  thereby  converting  pressure  energy  into  kinetic 
energy.  For  the  two-phase  nozzles  being  investigated  in  this  program,  the 
fluid  is  a  two-phase  mixture  consisting  of  hot  liquid  droplets  suspended  in  a 
gas.  Because  the  hot  liquid  acts  as  a  heat  source  in  a  two-phase  nozzle, 
part  of  the  thermal  energy  of  the  liquid  is  transferred  efficiently  to  the 
gas  phase,  and  the  resulting  two-phase  enthalpy  is  then  converted  into  kinetic 
energy. 

As  an  example  of  flow  in  a  two-phase  nozzle,  small  droplets  of  a  hot 
liquid  which  has  a  low  vapor  pressure  (such  as  mineral-oil  hydraulic  fluids) 
can  be  mixed  at  the  nozzle  inlet  with  a  liquid  that  has  a  higher  vapor 
pressure  (such  as  water).  If  the  hot  liquid  is  at  a  sufficiently  high 
temperature,  the  cold  liquid  with  its  high  vapor  pressure  can  be  vaporized  by 
the  heat  from  the  hot  liquid,  and  the  resulting  two-phase  mixture  is  then 
expanded  through  a  nozzle  to  low  pressure  and  high  velocity.  The  high-velocity 
vapor  drags  the  oil  droplets  along  to  a  high  velocity,  and  the  presence  of 
the  finely  dispersed  liquid  droplets  means  that  the  same  level  of  momentum 
flux  is  obtained  but  with  a  much  lower  gas  velocity  (or  spouting  velocity)  when 
compared  to  a  single-phase  gas  nozzle,  because  the  average  molecular  weight 
of  the  two-phase  mixture  is  higher  than  that  of  the  gas.  For  the  prime  mover 
application,  for  example,  two-phase  nozzles  can  be  integrated  with  a  two-phase 
turbine  for  efficient  operation  at  much  lower  rotor  speeds  than  a  turbine 
with  a  single-phase  gas  as  the  working  fluid.  Hence,  the  two-phase  turbine 
can  offer  high  torque  at  low  rotor  speed. 
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Flow  in  two-phase  nozzles  involves  changes  in  the  gas  velocity  and 
temperature  that  are  brought  about  by  the  gas-droplet  interactions  of 
viscous  drag  and  heat  transfer.  Because  of  the  short  residence  time  of  the 
two-phase  mixture  in  the  nozzle,  considerable  departures  from  velocity  and 
temperature  equilibrium  between  the  two  phases  can  occur  in  the  nozzle. 

Hence,  nonequilibrium  effects  should  be  included  in  analytical  models  of 
two-phase  nozzle  flows. 

Previous  work  on  two-phase  nozzles  can  be  divided  into  three  general 
areas:  1)  gas-particle,  2)  gas-liquid,  and  3)  vapor-liquid  flows.  Gas-particle 
nozzles  are  applied  mainly  to  metallized  solid  propellant  flows  in  rocket 
motors  and  are  the  simplest  to  analyze.  Gas-liquid  nozzles  are  more  complex, 
primarily  because  of  droplet  effects,  but  are  of  interest  in  a  wide  variety  of 
applications,  including  low-speed  power  turbines,  gas-driven  jet  pumps,  MHD 
power  generators,  and  atomization  and  spraying  equipment.  Vapor-liquid  nozzle 
flows  are  the  most  complex  because  of  the  interphase  transfer  mechanisms  and 
find  application  in  such  areas  as  geothermal  power  turbines  and  high-quality 
steam-water  flows  in  steam  turbines. 

1.1.1  Physical  Models  and  Solutions 

A  number  of  previous  investigators  have  proposed  physical  models  for 
two-phase  flows  in  nozzles,  and  several  of  these  models  are  discussed  in  the 
following  section.  Table  1.1  summarizes  the  important  aspects  of  these  models, 
all  of  which  are  one -dimensional  (except  for  Crowe's  model).  Generally,  the 
early  models  neglected  droplet  breakup  for  one  of  two  reasons:  either  to 
simplify  the  analysis,  or  to  reflect  the  fact  that  the  analysis  was  applicable 
to  constant-diameter  solid  spheres  in  rocket  exhaust  nozzles.  Most  of  the 
models  listed  also  neglect  shock  effects,  in  part  either  to  simplify  the 
analysis  or  to  reflect  the  fact  that  shock  effects  are  important  only  at  drop 
sizes  significantly  smaller  than  those  anticipated  for  their  application. 

As  seen  in  Table  1.1,  Tangren  et  al.  (Ref.  1)  conducted  one  of  the  early 
investigations  of  two-phase  flows  in  two-phase  nozzles,  and  nozzle  character¬ 
istics  for  isothermal  frictionless  flow  were  calculated.  The  isentropic  flow 
was  also  assumed  to  be  homogeneous,  which  assumption  means  that  differences 
in  velocity  and  temperature  were  neglected.  The  qualitative  behavior  of  the 
homogeneous  nozzle  flow  and  the  methods  of  solution  for  the  homogeneous  model 
are  very  similar  to  the  classic  closed-form  analytical  solution  for  the 
single-phase  nozzle  discussed  by  Shapiro  (Ref.  2).  The  Tangren  analysis 
applies  to  bubbly  flow,  but  should  also  be  valid  for  dispersed  (or  spray) 
flow  because  almost  all  of  the  mass  flux  is  also  due  to  the  liquid.  Further, 
when  the  liquid  comprises  most  of  the  mass  in  the  flow,  the  liquid  will  act 
as  a  heat  source  for  Che  expanding  gas,  and  the  expanding  two-phase  flow  may 
be  nearly  isothermal.  This  depends  on  the  droplet  size  which  affects  the 
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heat  transfer  rate  in  terms  of  both  the  heat  transfer  area  and  the  heat 
transfer  coefficient  between  the  gas  and  the  droplets.  The  Tangren  analysis 
produces  a  local  Mach  number  which  is  a  function  of  the  ratio  of  the  local 
pressure  to  the  reservoir  (or  inlet)  pressure,  and  the  throat  area  which  will 
pass  a  given  amount  of  the  two-phase  mixture.  However,  the  homogenous  theory 
does  not  account  for  differences  in  velocity  and  temperature  between  the  two 
phases,  and  because  differences  which  occur  in  the  real  nozzle  flow  lead  to 
nozzle  inefficiencies,  a  prediction  based  on  the  homogeneous  model  would  be 
optimistic. 

As  indicated  in  Table  1.1,  work  by  Kliegel  (Ref.  3)  took  into  consideration 
differences  in  velocity  and  temperature  between  the  two  phases.  The  analysis 
was  applied  to  gas-solid  flows  in  rocket  exhaust  nozzles,  where  the  size  of  the 
solid  particle  is  both  constant  and  known.  In  this  case,  too,  a  closed-form 
analytical  solution  was  obtained.  Kliegel  assumed  that  the  nozzle  shape  was 
such  that  the  ratio  of  particle  (or  droplet)  velocity  to  the  local  gas  velocity 
was  constant  and  that  the  lag  of  the  particle  temperature  is  directly  proportional 
to  the  lag  of  the  gas  temperature.  (The  lag  of  the  temperature  is  simply  the 
temperature  relative  to  the  reservoir  conditions  of  the  flow.)  These  nozzles 
are  known  as  "constant-fractional-lag  nozzles".  If  the  volume  fraction  of  the 
particles  is  neglected,  then  the  flow  equations  for  these  nozzles  reduce  to  the 
form  that  applies  to  the  flow  of  a  perfect  gas  which  has  an  apparent  specific 
heat  and  a  modified  Mach  number.  The  shape  of  the  nozzle  which  results  for  the 
assumption  of  constant  fractional  lag  is  arbitrary,  but  no  more  arbitrary  than 
other  shapes.  And  the  simplifications  thereby  possible  in  the  analysis  allow 
insight  into  the  behavior  of  the  two-phase  flow. 

Rudinger  (Ref.  4)  extended  Kliegel 's  concept  of  the  constant-fractional- 
lag  nozzle  to  high  particle  loadings  where  the  particle  volume  must  be 
considered.  For  this  situation,  computerized  solutions  are  necessary  because 
of  the  complicated  analysis.  Under  certain  flow  conditions,  the  computerized 
solutions  indicate,  however,  that  the  gas  and  particle  (or  droplet)  tempera¬ 
tures  remain  essentially  constant  throughout  the  length  of  the  nozzle  during 
certain  operating  conditions',  namely  those  with  high  loading  ratios  (which  is 
the  ratio  of  the  liquid  droplet/particle  mass  flow  rate  to  the  gas  flow 
rate).  These  results  are  plausible  because  isothermal  flow  implies  an 
isentropic  expansion  exponent  of  unity,  a  value  which  is  rapidly  approached 
as  the  loading  ratio  increases.  For  those  cases  where  the  flow  can  be 
assumed  to  be  isothermal,  then  the  equations  result  in  a  closed-form  analytical 
solution  convenient  for  hand  calculation  of  velocity  and  flow  area  as  a 
function  of  pressure  along  the  nozzle. 

Netzer  (Ref.  5)  included  the  effects  of  droplet  breakup  by  using  the  Weber 
number  breakup  criterion,  although  in  his  analysis,  friction  was  neglected. 

For  short  nozzles,  friction  effects  are  negligible;  however,  in  long  nozzles, 
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friction  effects  are  significant,  and  velocities  were  measured  that  were  as 
much  as  25  percent  below  those  predicted  by  this  frictionless  model.  Netzer's 
work  is  further  limited  by  the  fact  that  it  allows  for  no  temperature  difference 
between  the  two  phases. 

Hultberg  and  Soo  (Ref.  6)  investigated  two-phase  mixtures  of  gas  and 
solid  particles  (of  constant  diameter).  Their  model  included  velocity 
slip  and  heat  transfer  between  the  particles  and  the  gas  phase.  Iterative 
numerical  methods  were  used  to  solve  the  governing  equations  using  a  prescribed 
nozzle  profile,  and  the  convergence  of  the  solution  was  checked.  Because  the 
model  was  formulated  for  the  flow  of  a  metallized  propellant  containing 
alumina  or  magnesia  (solid)  particles  in  a  rocket  nozzle,  the  effects  of 
droplet  breakup  were  not  considered. 

Elliott  and  Weinberg  (Refs.  7  and  8)  synthesized  a  complex  computerized 
model  for  two-phase  nozzles  which  allowed  for  temperature,  as  well  as  velocity, 
differences  between  the  two  phases  and  included  droplet  breakup  and  friction 
effects  for  long  nozzles.  The  only  major  effect  not  considered  in  their  analysis 
was  that  of  shock.  This  model  was  used  (in  Ref.  8)  to  study  nozzle  designs  for 
an  MHD  application  in  which  predictions  were  compared  with  experimental  nitrogen- 
water  results.  For  this  application,  the  measured  bulk  velocities  and  flow 
rates  agreed  with  those  predicted  by  the  analysis  within  5  percent. 

Recent  work  by  Crowe  et  al .  (Ref.  9)  for  spray  cooling  resulted  in  a 
two-dimensional  model  which  includes  coupling  of  mass,  momentum,  and  energy 
between  the  two  phases,  although  droplet  breakup  was  not  included.  The  flow 
field  is  divided  into  a  series  of  cells,  and  finite-difference  equations  for 
mass,  momentum,  and  energy  conservation  are  written  for  the  gas  phase,  while 
the  droplet  effects  are  treated  as  source  terms  in  the  gas  flow  equations.  As 
a  numerical  example,  Crowe  analyzed  the  case  of  a  spray  issuing  into  a  moving 
gas  stream,  and  obtained  predictions  for  the  temperature  and  velocity  flow 
fields  of  the  gas  and  droplet  phases. 

Comfort  et  al .  (Ref.  10)  developed  a  one-dimensional  analytical  model  of 
two-phase  nozzle  flows  which  can  be  applied  to  steam-water  (one-component) 
mixtures  from  geothermal  heat  sources.  With  drop  diameters  of  less  than  10 
microns  (0.0004  in),  shock  effects  in  the  flow  can  be  important,  and  accordingly, 
those  were  included  in  the  analysis.  The  importance  of  accounting  for  two-phase 
shock  phenomena  increases  as  the  droplet  size  becomes  smaller,  because  the 
mixture  begins  to  behave  as  a  homogeneous  dense  gas.  Design  considerations  for 
the  geothermal  application  generally  require  that  the  velocities  in  such 
devices  as  two-phase  turbines  be  supersonic  relative  to  both  the  two-phase 
choking  velocity  and  the  vapor  choking  velocity.  In  the  Comfort  analysis, 
thermal  equilibrium  between  the  two  phases  was  assumed,  because  the  heat 
transfer  rate  is  very  high  when  the  drops  are  very  small.  For  drop  sizes  which 
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may  be  encountered  in  a  two-phase  nozzle  for  the  marine  propulsion  applications, 
however,  thermal  equilibrium  may  not  exist,  and  shock  effects  will  be  less 
significant.  Numerical  analysis  is  required  to  solve  Comfort's  model  which  is 
capable  of  predicting  pressure,  velocity,  droplet  diameter,  loading  ratio,  and 
normal  shock  waves  as  a  function  of  axial  distance  through  the  nozzle. 

1.1.2  Experimental  Data 

A  number  of  investigators  have  conducted  experiments  with  two-phase  nozzles, 
and  selected  parts  of  these  efforts  are  summarized  in  this  section.  Elliott 
and  Weinberg  (Ref.  8)  summarized  previous  work  as  of  1968  and  also  conducted 
experiments  with  room-temperature  nitrogen  and  water  mixtures  in  a  large-scale 
1.32  m  (52  inch)  long  nozzle  for  an  MHD  application.  Their  inlet  pressures 
ranged  from  10.3  to  15.9  bar  (150  to  230  psi),  and  their  expansion  pressure  ratio 
ranged  from  10:1  to  16:1.  With  a  loading  ratio  (liquid/gas)  of  40,  the  predicted 
temperature  drop  of  the  liquid  was  approximately  0.8  C  (1.5  F)  while  that  of  the 
gas  was  6  C  (11  F).  Without  the  water,  the  isentropic  expansion  of  the  nitrogen 
would  have  resulted  in  a  temperature  drop  of  139  C  (250  F).  Hence,  in  the 
two-phase  expansion,  the  liquid  acted  as  a  heat  source  for  the  gas,  with  liquid 
enthalpy  being  converted  into  gas  enthalpy.  Experimentally,  axial  pressure 
profiles  and  bulk  exit  velocities  were  measured  and  compared  with  the  analytical 
predictions.  The  flow  rate  data  generally  were  in  the  range  between  that 
predicted  for  flow  with  friction  and  that  for  isentropic  flow  (which  was  8 
percent  lower).  The  isentropic  prediction  is  lower  because  the  gas  velocity  is 
the  lowest  in  the  absence  of  slip,  and  hence  the  gas  flow  rate  is  lower.  And 
with  the  loading  ratio  fixed,  then  the  total  flow  rate  is  lower.  Elliott  and 
Weinberg  attributed  the  slightly  higher  predictions  of  flow  rate  to  the  under¬ 
estimation  of  the  droplet  drag.  Solid  sphere  drag  data  was  used,  whereas  the 
actual  liquid  droplets  become  flatter  and  offer  more  aerodynamic  resistance 
than  do  spherical  droplets.  Exit  velocities  were  also  measured  and  are  about  3 
to  5  percent  higher  than  those  predicted  for  flow  with  friction  -  the  difference 
is  again  attributed  to  underestimation  of  droplet  drag.  The  isentropic  velocity 
is  13  to  18  percent  higher  than  that  for  the  prediction  of  flow  with  friction. 
Tests  were  also  run  in  the  same  nozzle  using  Freon  1301  (CBrF^)  as  a  liquid 
at  the  nozzle  inlet  in  place  of  the  nitrogen  to  verify  that  the  gas  phase  could 
be  produced  by  contact  vaporization  at  the  inlet  of  a  two-phase  nozzle. 

Because  the  data  agreed  as  well  with  the  theory  as  in  the  nitrogen-water  tests 
(where  the  nitrogen  was  already  a  gas  at  the  nozzle  inlet),  Elliott  and  Weinberg 
concluded  that  complete  vaporization  of  the  Freon  1301  did  occur. 

Elliott  also  reported  (Ref.  7)  the  results  of  earlier  water-air  experiments 
using  a  0.15-m  (6-in.)  long  converging-diverging  nozzle  at  inlet  conditions  of 
12  C  (54  F)  and  35.4  bar  (514  psia).  These  results  were  compared  with  those  of 
their  later  analytical  model  (Ref.  8).  The  measured  flow  rates  were  only  1  to 
2  percent  lower  than  those  predicted  for  flow  with  friction  and  the  measured 
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bulk  exit  velocities  straddled  by  +  3  percent  the  prediction  for  flow  with 
friction.  For  this  short  nozzle,  the  effect  of  friction  was  to  decrease  the 
flow  rate  prediction  based  on  the  frictionless  prediction  criterion  by  3 
percent,  and  the  bulk  exit  velocity  prediction  by  6  percent. 

Starkman  et  al .  (Ref.  11)  measured  flow  rates  of  one-component,  two-phase 
(water-steam)  mixtures  expanding  through  a  converging-diverging  deLaval  nozzle, 
with  the  inlet  quality  ranging  from  zero  to  20  percent  and  inlet  pressures  up 
to  69  bar  (1000  psia) .  The  data  were  compared  with  the  results  of  three  simple 
models:  isentropic  expansion  with  thermal  equilibrium  but  no  slip;  frozen 
composition  with  no  slip  or  heat  transfer;  and  slip  flow  with  limited  momentum 
exchange  between  the  phases.  The  measured  flow  rates  agreed  satisfactorily 
with  the  results  from  all  three  models,  except  at  low  inlet  qualities  (where 
the  mixture  is  mostly  saturated  liquid).  When  the  nozzles  were  run  at 
overexpanded  conditions,  it  was  observed  that  there  was  a  two-phase  "shock" 
which  was  much  weaker  and  less  distinct  than  that  found  in  single-phase  gas 
flows.  Schrock  et  al .  (Ref.  12)  extended  the  work  of  Starkman  et  al .  and 
examined  the  effect  of  inlet  subcooling  on  the  flashing  flow  of  water  in 
converging-diverging  nozzles.  A  two-step  model  based  on  a  nucleating  delay 
with  a  discontinuous  transition  to  a  frozen  composition  gave  reasonable 
predictions  of  flow  rates  and  pressure  profiles  in  the  convergent  section. 

For  the  geothermal  application  of  hot  water  deposits,  Alger  (Ref.  13) 
tested  a  small  nozzle  (0. l-m(4-inches)long)  with  mixtures  of  steam-water  of  19 
percent  inlet  quality  at  24.8  bar  (360  psia)  and  224  C  (435  F) .  This  nozzle 
could  be  incorporated  into  a  two-phase  impulse  turbine,  although,  to  achieve  a 
maximum  efficiency  turbine,  droplet  sizes  of  less  than  4  microns  (0.00016  in) 
would  be  required  so  that  the  droplets  will  follow  the  vapor  streamlines  and 
not  collide  with  the  blade  walls.  The  thrust  of  Alger's  work  was  to  measure 
droplet  phase  characteristics  (primarily  drop  sizes),  and  the  results  suggest 
that  the  droplets  break  up  at  diameters  which  are  less  than  about  one  third 
of  the  diameter  predicted  by  a  critical  Weber  number  of  6.  (The  Weber  number 
is  defined  as  the  dimensionless  ratio  of  the  aerodynamic  force  to  the  surface 
tension  force  on  the  droplets,  and  the  critical  Weber  number  is  the  value  at 
which  droplet  breakup  occurs.) 

Cerini  (Ref.  14)  tested  a  0.63-m  (25-inch)  long  steam-water  nozzle  for 
a  geothermal  application  where  inlet  conditions  were  5.5  bar  (80  psia)  and  5 
percent  quality.  The  measured  flow  rates  agreed  to  +3  percent  with  the 
predictions  from  the  computerized  model  with  friction  developed  by  Elliot  and 
Weinberg  (Ref.  8). 

Most  of  the  recent  experimental  work  on  two-phase  nozzles  has  concentrated 
on  geothermal  applications,  where  the  steam-water  mixtures  are  one-component 
and  mass-transfer  effects  are  very  important.  However,  two-phase  nozzles  with 
two-component  flows  are  of  interest  for  applications  in  marine  power  and  thrust 
augmentation  in  turbofan  engines. 
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1.1.3  Two-Phase  Nozzle  Applications 

Two-phase  nozzles  are  of  interest  in  a  number  of  fields,  and  Table  1.2 
lists  several  of  these  applications.  As  can  be  seen,  two-phase  nozzles  are 
being  used  in  such  areas  as  propulsion,  energy  conversion,  and  tranaportation. 
Two-phase  nozzles  are  used  in  two-phase  turbines  which,  because  of  the  higher 
density  of  the  two-phase  working  fluid,  have  lower  rotor  speeds  than  do  tur¬ 
bines  using  gas  as  the  working  fluid.  With  the  lower  rotor  speeds,  two-phase 
turbines  (Ref.  15)  offer  the  possibility  of  either  a  smaller  gearbox  or  no 
gearbox  in  the  marine  and  automotive  propulsion  applications.  With  a  two-phase 
nozzle,  the  liquid/gas  ratio  can  be  varied  to  control  the  nozzle  exit  velocity, 
a  feature  which  is  unique  to  two-phase  nozzles.  As  a  result,  the  turbine  speed 
can  be  varied  at  nearly  constant  efficiency.  At  part-load  conditions,  the 
efficiency  of  a  two-phase  turbine  cycle  can  be  higher  than  that  of  a  single¬ 
phase  steam  turbine,  because  the  efficiency  of  a  two-phase  turbine  cycle  is 
nearly  flat  over  a  large  portion  of  its  turbine  speed  range.  Further,  the 
usual  boiler  could  be  eliminated  because  a  single-phase  liquid  heat  exchanger 
could  be  used  to  heat  the  liquid  which,  in  turn,  would  be  mixed  with  a  high- 
vapor-pressure  liquid.  The  resulting  vapor  would  be  expanded  through  the 
the  nozzle  carrying  with  it  the  hot  liquid  drops  (of  the  low-vapor-pressure 
fluid) . 

Another  application  of  great  current  interest  is  geothermal  power, 
where  the  steam  is  frequently  unsaturated,  and  where  potential  erosion  problems 
play  a  decisive  role  in  the  design  of  the  turbine.  With  the  higher  density  of 
the  two-phase  mixture,  the  turbine  rotor  speed  is  lower  than  in  a  gas  turbine, 
and  erosion  problems  would  be  less  severe. 

Water  injection  in  the  fan  exhaust  nozzle  of  a  turbofan  engine  can  be  used 
to  augment  the  thrust.  Work  on  this  two-phase  nozzle  (Ref.  16)  using  a  water- 
augmented  turbofan  engine  for  high-speed  ships  showed  a  theoretical  augmentation 
in  thrust  of  380  percent  at  12.9  m/sec  (25  knots)  and  of  90  percent  at  51.4 
m/sec  (100  knots). 

Work  by  Elliot  and  Weinberg  (Ref.  8)  was  directed  toward  MHD  power  generation, 
where  liquid  lithium  and  cesium  vapor  are  the  two  components.  Ahmad  and  Hays 
(Ref.  18)  have  investigated  the  use  of  a  two-phase  turbine  system  to  recover 
waste  heat  in  a  bottoming  cycle  of  a  stationary  diesel  engine,  and  a  fuel 
saving  of  as  much  as  30  percent  is  calculated. 


1.2  Identification  of  Governing  Flow  Parameters 

Preparatory  to  formulating  a  working  theory  of  two-phase  nozzle  flows,  the 
governing  flow  parameters  and  their  ranges  of  values  were  identified.  Liquid/gas 
combinations  for  the  various  two-phase  nozzle  applications  were  examined,  and 
suitable  combinations  are  discussed  for  selected  applications. 
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1.2.1  Governing  Flow  Parameters 

The  principal  governing  flow  parameters  were  identified  on  the  basis  of  the 
previous  investigations  that  were  described  in  Section  1.1.  Fig.  1.1  is  a 
schematic  of  a  converging-diverging  nozzle  with  the  governing  flow  parameters: 
the  velocity(v),  mass  flow  rate(m),  and  temperature(T)  of  both  the  gas-phase 
(subscript  g)  and  liquid-phase  (subscript  i)  streams,  and  the  droplet  diameter(D). 
The  thermophysical  properties  of  each  stream  must  also  be  known.  For  the  gas 
phase,  these  properties  include  the  density  (Pg)>  specific  heat  (Cpg), 
thermal  conductivity  (k  ) ,  and  viscosity  (p^).  For  the  liquid  phase,  these 
properties  include  the  density  (p^),  specific  heat  (C),  and  the  surface 
tension  (o)  for  use  in  estimating  droplet  sizes.  The  local  differences  in 
velocity  and  temperature  between  the  two  phases  affect  the  interphase  transfer 
mechanisms  of  momentum  and  energy.  The  velocity  difference  affects  the 
droplet  drag,  the  heat  transfer  coefficient  between  the  liquid  droplet  and 
the  gas,  and  the  droplet  breakup.  The  temperature  difference  is  the  driving 
force  for  heat  transfer  between  the  liquid  and  the  gas.  The  flow  rates  of 
each  phase  determine  the  loading  ratio  (m^/m  )  of  the  nozzle  and 
the  mass  flux  of  each  phase  in  the  nozzle.  Tne  droplet  diameter  affects  the 
coefficients  of  drag  and  heat  transfer  and  is  obtained  from  the  Ueber  droplet 
breakup  criterion  which  in  turn  depends  on  the  velocity  difference  between 
the  two  phases. 

The  range  of  parameters  considered  in  previous  investigations  of  two-phase 
nozzles  is  shown  in  Table  1.3  in  terms  of  the  gas/liquid  combination,  the  design- 
point  temperature  and  pressure,  the  liquid  mass  and  volume  fractions,  and  the 
droplet  diameter  (when  indicated  by  the  investigator).  In  a  given  investigation, 
the  analytical  and  experimental  aspects  generally  covered  the  same  range  of 
conditions,  except  when  testing  was  done  with  fluids  that  are  easier  to  handle. 

For  example,  Elliot  and  Weinberg  (Ref.  8)  examined  the  application  of  MHD  power 
generation  with  liquid  metals  (lithium  as  the  liquid  and  cesium  as  the  gas 
phase)  at  inlet  conditions  of  982  C  (1800  F)  and  9.44  bar  (137  psia),  but  their 
testing  was  done  with  mixtures  of  nitrogen  and  water  at  inlet  conditions  of  16  C 
(60  F)  and  10.3  bar  (ISO  psia). 

In  related  investigations  at  the  Lawrence  Livermore  Laboratory,  Alger 
(Ref.  13)  and  Comfort  and  Beadle  (Ref.  17)  examined  the  geothermal  application, 
and  both  testing  and  analysis  were  performed  with  steam/water  mixtures  at  the 
actual  operating  conditions  [224  C  (435  F)  and  24.8  bar  (360  psia)]  of  a  geothermal 
heat  source.  The  liquid  mass  fraction  of  these  investigations  ranged  from 
approximately  0.6  to  1.0  (which  is  all  water). 
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The  expected  droplet  diameter  also  covered  a  wide  range,  from  l  to  10  microns 
(3. 9x10“'*  to  3.9x10“^  in.)  in  Alger's  investigation  to  large  droplets  of 
1270  microns  (0.050  in.)  in  Elliot  and  Weinberg's  MHD  investigation.  Analyti¬ 
cally,  Alger  showed  that  a  large  increase  in  the  efficiency  of  a  total-flow 
impulse  turbine  occurs  when  the  droplet  size  is  decreased  because  if  the 
droplets  are  too  large,  they  do  not  follow  the  vapor  streamlines  within  the 
turbine  blade  passages  and  collide  with  the  turbine  walls,  thereby  resulting  in 
reduced  momentum  exchange  and  liquid  film  pumping  losses.  However,  in  the 
separating  type  of  two-phase  turbine,  smaller  droplets  can  adversely  affect  the 
separation  efficiency,  because  these  smaller  droplets  are  more  difficult  to 
separate  from  the  gas  in  a  centrifugal  field. 

Analytically,  Elliot  and  Hays  (Ref.  15)  reviewed  two-phase  nozzles  and 
turbines  for  the  power  turbine  application,  such  as  marine  propulsion.  Here, 
the  nozzle  inlet  conditions  for  the  mixture  of  steam  and  polyphenyl  ether 
(which  is  a  vacuum  pump  oil)  were  327  C  (620  F)  and  100  bar  (1450  psia),  and 
the  expected  drop  size  approximately  25  microns  (9.8xl0~^  in.).  In  a  related 
investigation  (Ref.  18)  for  a  bottoming  cycle,  the  design  conditions  at  the 
nozzle  inlet  were  349C  (660  F)  and  163  bar  (2365  psia)  for  mixture  of  steam 
and  Dowtherm  A.  As  a  point  of  reference,  the  critical  point  of  steam  is  392  C 
(705  F),  above  which  point  there  is  no  phase  change  from  a  liquid  to  a  vapor 
and  the  steam  would  act  increasingly  as  a  heavy  liquid,  thereby  losing  the 
acceleration  effects  of  the  compressible  gas. 

The  ranges  of  the  governing  flow  parameters  for  a  two-phase  nozzle  to  be 
used  in  a  two-phase  turbine  can  be  estimated  as  follows.  If  a  mixture  of  steam 
and  heat  transfer  oil  is  selected  as  the  working  fluid,  then  the  siaximum 
working  temperature  will  be  about  370  C  (700  F),  because  the  heat  transfer  oil 
will  begin  to  decompose  at  higher  temperatures  according  to  the  manufacturer's 
data,  and  because  the  critical  temperature  of  steam  is  392  C  (705  F).  Hence, 
the  inlet  temperature  will  be  in  the  vicinity  of  316  C  (600  F).  As  to  the 
exit  temperature,  because  of  the  anticipated  high  loading  ratios,  the  equili¬ 
brium  exit  temperature  of  the  mixture  will  be  only  slightly  lower  than  the 
inlet  temperature.  The  homogeneous  flow  model  can  be  used  to  show  this,  as 
indicated  by  the  results  in  Fig.  1.2  for  typical  mixtures  of  steam  and  Dowtherm 
A  at  316  C  (600  F)  and  two  pressure  levels,  96.6  and  69  bars  (1400  and  1000 
psia).  Here,  the  effect  of  loading  ratio  on  the  temperature  drop  of  a  homogene¬ 
ous  two-phase  mixture  expanding  through  a  converging-diverging  nozzle  is  shown, 
at  two  typical  values  of  pressure  ratio.  For  loading  ratios  larger  than  10, 
the  temperature  drop  of  the  homogeneous  mixture  is  less  than  3  percent  at  an 
expansion  ratio  of  2,  and  less  than  9  percent  at  a  ratio  of  10.  Hence,  the 
homogeneous  temperature  will  be  nearly  isothermal  at  some  of  the  conditions 
of  interest. 
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Zn  Actuality,  however,  heat  tranafer  ia  a  rate  process,  and  the  droplet 
size  will  alao  have  a  large  effect  on  whether  the  expansion  ia  nearly  iaothermal 
becauae  the  droplet  aize  will  affect  both  the  value  of  the  heat  transfer 
coefficient  between  the  gaa  and  the  droplet  and  the  amount  of  heat  tranafer 
area  of  the  droplet,  aa  will  be  aeen  later.  Fig.  1.3  ia  a  Temperature-Entropy 
diagram  of  the  cycle  paths  for  one-phase  and  two-phase  turbines.  Aa  indicated 
by  the  cycle  path  above  the  cross-hatched  area,  the  use  of  a  two-phase  nozzle 
can  result  in  an  expanaion  that  is  nearly  isothermal,  because  the  liquid  acts 
as  a  heat  source  for  the  expanding  gas  (which  otherwise  would  be  cooled  in  a 
single-phase  expansion).  This  heat  transfer  produces  additional  work  which  can 
be  extracted  by  the  two-phase  turbine  as  represented  by  the  heat  in  the  integral 
of  TdS  shown  as  the  cross-hatched  area  in  Fig.  1.3.  The  saturation  pressure  of 
steam  at  316  C  (600  F)  is  107  bar  (1550  pai),  and  hence,  the  nozzle  inlet 
pressure  should  be  established  at  this  level  so  the  area  enclosed  by  the 
Temperature-Entropy  diagram  (as  shown  schematically  in  Fig.  1.3  for  the 
entire  cycle)  will  be  large.  The  exit  pressure  should  be  as  low  as  possible  s.' 
that  the  conversion  of  enthalpy  to  kinetic  energy  will  be  as  large  as  possible. 
The  condensing  temperature  will  be  typically  43  C  (100  F),  where  the  saturation 
pressure  for  steam  is  0.088  bar  (1.27  psia).  Hence  the  nozzle  exit  pressure 
must  be  higher  than  that  to  allow  for  pressure  drops  through  the  system,  and 
will  be  in  the  area  of  1  bar  (14.5  psia),  which  is  low  enough  for  nozzle 
expansion  and  yet  high  enough  to  allow  for  pressure  in  the  flow  path  to  the 
condenser. 

Based  on  work  by  Elliot  and  Hays  (Ref.  15),  the  loading  ratio  (mi/mg) 
of  the  nozzle  will  range  from  about  10  to  50.  At  a  loading  ratio  of  50,  the 
speed  of  the  two-phase  turbine  will  be  about  one-third  that  of  the  all-gas 
turbine,  with  only  a  slight  drop  in  cycle  efficiency.  Larger  values  of  the 
loading  ratio  result  in  a  significant  drop  in  cycle  efficiency.  At  the  other 
end,  a  loading  ratio  of  10  is  near  the  lower  limit,  because  ratios  lower  than 
that  result  in  turbine  speeds  that  approach  the  all-gas  turbine  speed,  and 
the  low-speed  advantage  of  a  two-phase  turbine  would  then  be  lost. 

The  expected  droplet  size  can  be  deduced  from  the  Weber  number  droplet 
breakup  criterion.  Fig.  1.4  shows  calculated  results  for  the  diameter  of 
droplets  of  Dow therm  A  in  steam  as  a  function  of  steam  velocity,  idien  the 
maximum  Weber  number  is  equal  to  6.  (The  Weber  number  is  a  dimensionless 
criterion  which  is  the  ratio  of  the  aerodynamic  pressure  forces  which  are 
trying  to  deform  and  fragment  the  droplet  and  the  surface  tension  forces  which 
are  trying  to  maintain  the  droplet  in  spherical  form.  The  critical  value  of 
the  Weber  number  is  experimentally  determined.)  At  the  nozzle  inlet  where  the 
gas  (steam)  velocity  will  be  approximately  30.5  m/sec  (100  ft/sec)  and  droplet 
velocities  will  be  in  the  range  of  18.3  to  24.4  m/sec  (60  to  80  ft/sec),  the 
expected  drop  size  will  be  in  the  range  of  50  to  180  microns  (0.002  to  0.007  in.) 
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1.2.2  Liquid/Gas  Combinations 

A  liquid/gas  combination  suitable  for  each  two-phase  nozzle  application 
must  be  selected.  In  general,  there  are  two  types  of  two-phase  liquid/gas 
flows.  If  the  gas  phase  is  the  vapor  of  the  liquid  phase,  then  the  flow  is 
"one-component"  flow,  and  an  example  of  this  type  is  steam  and  water  flow.  If 
the  gas  phase  is  of  a  different  chemical  species  than  that  of  the  liquid,  then 
the  flow  is  "two-component"  flow,  and  an  example  of  this  would  be  a  mixture  of 
steam  and  an  organic  liquid  with  a  low  vapor  pressure.  The  expansion  of  a 
one-component  system  occurs  along  the  saturation  line,  while  the  expansion  of  a 
two-component,  system  can  be  nearly  isothermal  if  the  loading  ratio  is  high 
enough.  A  nearly  isothermal  expansion  is  advantageous  for  a  two-phase  nozzle 
in  a  two-phase  turbine.  For  a  two-component  system,  the  main  factors  in  the 
selection  of  the  liquid  are  sufficiently  high  density,  low  vapor  pressure,  and 
inertness,  while  the  main  requirements  for  the  gas  are  iimniscibility  with  the 
liquid  and  a  reasonable  condensing  pressure.  This  pressure  must  not  be  so  low 
as  to  require  a  huge  nozzle  exit  area,  and  not  be  so  high  as  to  require  high 
nozzle  exit  area. 

For  many  applications,  the  liquid/gas  combination  is  an  inherent  characteristic 
of  that  application:  the  steam/brine  mixture  for  the  geothermal  application 
because  the  thermal  energy  is  already  contained  in  that  mixture,  or  the  air/water 
mixture  for  thrust  augmentation  in  marine  propulsion  because  air  is  the  working 
fluid  of  the  gas  turbine  and  water  is  both  suitable  and  inexpensive  for  this 
open-cycle  process.  The  MHD  power  generation  application  requires  high  working 
temperatures  in  the  vicinity  of  982  C  (1800  F) ,  which  means  that  all  working 
fluids  except  for  the  liquid  metals  are  ruled  out.  Here,  the  liquid/gas 
combination  of  greatest  interest  (Ref.  8)  is  the  two-component  system  of  liquid 
lithium  and  cesium  vapor,  but  a  one-component  system  of  potassium  and 
potassium  vapor  is  also  of  interest  for  MHD.  Another  application  of  two-phase 
nozzles  (Ref.  7)  is  that  of  jet  pumps  (with  no  moving  parts)  for  rocket 
engines.  Here,  a  non-condensible  gas  (such  as  hydrazine  decomposition  products) 
is  used  to  pump  a  liquid  propellant  (such  as  nitric  acid,  hydrogen  peroxide, 

JP-5  jet  fuel,  or  hydrazine,  etc). 

Recently,  Ahmad  and  Hays  (Ref.  18)  have  studied  the  use  of  a  two-phase  turbine 
for  a  diesel  bottoming  cycle.  The  nozzle  inlet  temperature  is  349  C  (660  F), 
and  after  an  extensive  cycle  study,  a  mixture  of  Dowtherm  A  and  steam  was 
chosen  as  the  liquid/gas  combination.  The  use  of  a  liquid  with  a  lower  specific 
heat  than  Dowtherm  A  would  mean  a  higher  loading  ratio  and  hence  a  lower  nozzle 
exit  velocity,  but  the  fluid  with  the  lowest  specific  heat  that  they  found  was 
a  "Therminol"  fluid,  which  may  have  toxicity  problems.  The  thermophysical 
properties  (vapor  pressure,  specific  heat,  density,  viscosity,  and  thermal 
conductivity)  are  tabulated  as  a  function  of  temperature  in  Ref.  16  for  11 
fluids  of  potential  interest  for  an  application  at  the  316  C  (600  F)  level. 

Included  among  these  fluids  are  3  kinds  of  Therminol,  2  kinds  of  Dowtherm,  and 
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2  kinds  of  DowCorning  fluids.  (Another  fluid  of  possible  interest  is  polyphenyl 
ether,  and  three  chemical  firms  were  contacted:  Union  Carbide  does  not  make 
it,  Dow  Chemical  stopped  making  it  20  years  ago,  but  Monsanto  does  make  it  for 
use  in  vacuum  diffusion  pumps  and  sells  it  for  approximately  $1000/gal.)  These 
same  eleven  fluids  are  also  of  interest  in  the  marine  propulsion  application 
because  the  range  of  operating  temperatures  is  similar. 

Two-phase  turbines  are  being  considered  (Ref.  15)  for  marine  and  automotive 
propulsion  applications  because  of  their  low-speed  characteristics  which  result 
in  reduced  gearbox  requirements.  High-temperature  synthetic,  organic  heat- 
transfer  oils  (such  as  the  Dowtherm  family)  seem  to  offer  the  most  advantages 
as  the  liquid  for  this  application  where  steam  would  be  the  gas  phase.  These 
oils  have  a  density  only  about  15  percent  less  than  that  of  water  and  a  very 
low  vapor  pressure.  Their  specific  heat  is  in  the  range  of  0.5  cal/gm-C  (0.5 
BTU/lbm-F)  which  means  that  the  liquid  flow  rate  (and  hence  the  loading  ratio) 
will  be  high  enough  to  result  in  a  reasonably  low  nozzle  exit  (or  spouting) 
velocity.  Further,  most  of  these  oils  are  inert  and  have  a  low  level  of 
toxicity.  For  the  gas  phase,  water  is  a  well-known  working  fluid  which  has  a 
reasonable  vapor  pressure  characteristic,  and  is  immiscible  with  the  heat 
transfer  oils  under  consideration  for  the  liquid  phase.  In  general,  fluorocar¬ 
bons  are  not  suitable  for  use  as  the  gas  phase  for  three  reasons.  First,  as 
seen  in  Fig.  1-5,  their  vapor  pressures  are  significantly  higher  than  that  of 
water  at  the  same  pressure.  Second,  as  also  seen  in  Fig.  1.5,  their  critical 
temperatures  are  significantly  lower  than  that  of  water,  which  means  that  if 
operating  above  the  critical  temperature,  compressibility  effects  in  accelerat¬ 
ing  the  mixture  would  be  lost  because  the  fluorocarbon  would  behave  essentially 
as  a  heavy  liquid.  Third,  these  fluorocarbons  decompose  (into  toxic  substances) 
at  temperatures  above  177  C  (350  F),  and  this  decomposition  is  only  aggravated 
by  the  presence  of  foreign  substances  such  as  water  or  oils. 


1.3  Classification  of  Flow  Domains  With  Available  Theory 

The  flow  domains  considered  and  the  theories  used  in  previous  investigations 
of  two-phase  nozzles  were  examined  for  their  usefulness  in  classifying  and 
correlating  flow  domains.  Single-phase  (or  homogeneous)  theory  was  reviewed 
from  the  viewpoint  of  energy  conversion  in  a  nozzle  and  was  compared  to  models 
with  separate  phases  in  terms  of  their  applicability  to  these  flow  domains. 

Gas  dynamic  aspects  of  two-phase  nozzle  flows  were  examined  in  terms  of  the 
speed  of  sound  in  a  two-phase  mixture  and  of  choked  flow  at  the  nozzle  throat. 

1.3.1  Single-Phase  Nozzle  Theory 

For  some  operating  conditions  of  two-phase  nozzles,  the  flow  regime  will  be 
dispersed  flow,  where  the  liquid  phase  is  in  the  form  of  uniformly  dispersed, 
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small-diameter  droplets  suspended  in  the  gas.  In  many  ways,  the  behavior  of 
droplets  suspended  in  a  gas  is  similar  to  the  behavior  of  bubbly  flow  where  the 
small  bubbles  are  suspended  in  the  liquid.  In  bubbly  flow,  most  of  the  inertia 
is  in  the  liquid  phase,  and  the  drag  forces  on  the  gas  bubbles  are  large 
relative  to  the  momentum  forces.  As  a  result,  the  bubbles  follow  the  forced 
convection  motion  of  the  surrounding  liquid  very  closely.  For  dispersed  flow 
where  the  situation  is  reversed  with  the  liquid  droplets  suspended  in  the  gas, 
the  momentum  forces  are  much  larger  than  the  drag  forces,  and  the  droplets  do 
not  adjust  to  the  flow  of  the  surrounding  gas  as  fast.  For  this  reason, 
homogeneous  flow  models  are  not  expected  to  describe  the  situation  that  obtains 
in  dispersed  flow  as  well  as  in  bubbly  flow,  but  useful  understanding  of  the 
basic  flow  processes  can  still  be  obtained  from  a  homogeneous  model. 


The  theory  for  single-phase  incompressible  flow  is  valid  for  two  cases.  The 
first  case  is  when  the  nozzle  exit  pressure  is  only  slightly  less  than  the 
inlet  pressure  and  the  flow  (either  gas  or  liquid)  can  then  be  treated  as 
incompressible  because  density  effects  are  negligible.  The  second  case  is  when 
the  flow  is  all  liquid  (which  corresponds  to  a  loading  ratio  of  infinity),  and 
there  are  no  compressibility  effects.  Essentially,  the  continuity  equation  is 
then  used  to  calculate  nozzle  performance,  as  there  is  no  pressure  effect  and 
only  a  small  temperature  effect  on  liquid  density.  If  the  flow  were  mostly 
liquid  with  only  a  few  small  bubbles,  then  the  flow  would  be  expected  to  behave 
as  a  liquid  with  only  small  compressibility  effects. 


At  the  other  extreme  where  the  single-plase  flow  is  all  gas  (and  the  loading 
ratio  is  zero),  the  classical  isentropic  compressible  gas  theory  (Ref.  2), 
describes  the  flow  situation.  For  the  flow  situations  between  the  all-liquid 
and  the  all-gas  flows,  single-phase  isentropic  theory  has  been  developed  to 
describe  the  liquid/gas  mixture  in  terms  of  average  properties,  assuming 
thermal  equilibrium  (T^  *  Tg)  and  no  velocity  slip  (Vfc  ■  Vg).  The  main 
influences  on  the  behavior  of  such  a  mixture  are  the  compressibility  of  the  gas 
phase  and  the  inertia  of  the  liquid.  This  theory  (called  the  homogeneous  flow 
theory  and  given  in  Ref.  23,  for  example)  treats  the  liquid/gas  mixture  as  a 
pseudo  gas  which  is  then  shown  to  obey  the  well-known  one-dimensional  flow 
relationships  of  classical  gas  dynamics.  The  isentropic  expansion  exponent  of 
the  mixture,  r,  is: 


r 


Hr& 

1  ♦  V5 


(1.1) 


where  ig  is  the  isentropic  expansion  exponent  of  the  pure  gas,  r  is  the 
loading  ratio  (m^/nig),  and  6  is  the  ratio  of  the  specific  heat  of  the 
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liquid  to  the  specific  heat  of  the  gas  at  constant  pressure  (C/Cpg). 

For  the  all-liquid  extreme,  the  loading  ratio  (n)  is  infinite  and  r  *  1.  At 
the  all-gas  extreme,  the  loading  ratio  is  zero  and  r  *  y  .  Then  the  temperature 
and  density  after  an  isentropic  expansion  are  given  by  the  expressions: 


(1.2) 


where  the  subscript  zero  refers  to  nozzle  reservoir  (or  inlet)  conditions.  This 
expression  for  temperature  in  an  isentropic  expansion  was  used  to  obtain  the 
results  shown  earlier  in  Fig.  1-2.  The  nozzle  exit  velocity,  V2,  is  given  by 
the  expression: 


(1.3) 


where  Rg  is  the  gas  constant  of  the  pure  gas.  An  expression  for  the 
momentum  flux  at  the  nozzle  exit  can  be  obtained  using  continuity  and  the  above 
expressions  for  the  nozzle  exit  velocity  and  density: 

1  r-i ,  i/r 

f  -°>l  m°fe)!>-(fc)r  l(y 


Results  calculated  from  these  expressions  for  exit  velocity  and  momentum 
flux  for  mixtures  of  steam  and  Dowtherm  A  are  shown  in  Fig.  1.6  as  a  function 
of  the  loading  ratio.  Several  observations  can  be  made  from  these  results. 
First,  the  all-gas  flow  corresponds  to  a  loading  ratio  of  zero,  while  the 
all-liquid  extreme  is  approached  as  r  ♦  ®.  Second,  the  momentum  flux  is 
nearly  independent  of  loading  ratio,  because  the  amount  of  enthalpy  converted 
into  kinetic  energy  (in  the  form  of  momentum)  is  essentially  fixed  by  the 
expansion  ratio.  Third,  as  the  loading  ratio  increases,  the  nozzle  exit 
velocity  decreases.  Again,  this  iB  expected  because  the  momentum  flux  is 
nearly  constant,  and  the  density  of  the  homogeneous  two-phase  mixture  is 
increasi. ^  as  the  loading  ratio  is  increased. 

1.3.2  Comparison  of  Homogeneous  Model  With  Separate-Phase  Models 

A  flow  model  with  partially  separate  phases  was  developed  by  Rudinger 
(Ref.  4)  when  he  extended  Kliegel's  constant-fractional-lag  model  (Ref.  3)  to 
include  the  effect  of  the  particle  volume  fraction.  Using  the  equations  of 
continuity,  momentum,  energy,  and  the  heat  balance  on  a  particle  (or  droplet), 
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a  model  was  developed  and  computerized  numerical  solutions  were  obtained.  But 
these  solutions  indicated  that  the  gas  and  particle  temperatures  remained 
substantially  constant  throughout  the  nozzle  over  a  wide  range  of  operating 
conditions  of  interest.  With  the  assumption  of  isothermal  flow  (T^  *  Tg  m  TQ), 
then  the  energy  and  the  heat  balance  equations  are  eliminated,  and  a  particularly 
simple  result  is  obtained  for  the  gas  velocity  in  the  nozzle  as  a  function  of 
pressure: 


+  2RgTa  [W&L  ££_ 

1+rK  l  \P  /  l-e0 


(1.5) 


where  K  is  the  velocity  lag  (K  =  Va/Vg),  and  e  is  the  particle  volume 
fraction.  The  subscript  zero  refers  to  inlet  or  reservoir  conditions.  For  the 
special  case  when  no  particles  (or  droplets)  are  present  in  the  flow,  both  the 
loading  ratio  (r)  and  the  inlet  particle  volume  fraction  (eQ)  are  zero,  and 
the  isothermal  solution  of  Eq.  1.5  reduces  to: 

Vg2  “  2¥o  In  Pq  (1. 5a) 

P 

if  the  inlet  velocity  is  zero.  For  steam  at  260  C  (500  F),  the  calculated  exit 
velocity  from  Eq.  I. 5a  for  an  isothermal  expansion  is  502  m/sec  (1646  ft/sec) 
at  an  expansion  ratio  (PQ/P)  of  1.66,  and  1066  m/sec  (3496  ft/sec)  when  PQ/P  is 
10.  In  contrast,  for  the  isentropic  expansion  of  pure  steam  (in  the  absence  of 
liquid  oil  droplets),  the  homogeneous  model  of  Eq.  1.3  yields  an  exit  velocity 
of  488  m/sec  (1599  ft/sec)  at  an  expansion  ratio  of  1.66,  and  939  m/sec  (3079 
ft/sec)  when  PQ/P  is  10.  Hence,  the  isothermal  model  for  pure  steam  predicts 
nozzle  exit  velocities  that  are  5  to  10  percent  higher  than  the  homogeneous 
model.  Continuity  considerations  would  lead  one  to  expect  a  higher  velocity 
for  the  isothermal  case  where  the  temperature  is  higher  (and  hence  the  gas 
density  is  lower  at  the  same  pressure).  For  the  isothermal  model,  expressions 
for  the  other  variables  are: 

-L.  I_ 

1’e  1-eo  po 

and 

A_  .  Vg)oPo  /  l^j 
A0  VgP  '  1-e  f 

By  using  the  continuity  equations  of  the  gas  and  of  the  liquid,  the  following 
relationship  between  the  local  values  of  the  flow  variables  can  be  obtained: 


(1.6) 


(1.7) 
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The  numerical  relationship  between  the  variables  of  Eq.  1.8  is  shown  in  Fig.  1.7. 
For  example,  mixtures  of  saturated  steam  and  Dowtherm  A  at  260  C  (500  F)  and 
46.9  bar  (680  psia)  have  a  density  ratio  (P^/Pg)  of  approximately  36. 

Thus,  for  a  loading  ratio  of  20  (that  could  be  typical  for  a  marine  propulsion 
application),  the  particle  volume  fraction  is  over  0.40  for  the  slip  ratio  in  a 
typical  range  of  0.6  to  0.8.  In  other  words,  the  gas  volume  accounts  for  under 
40  percent  of  the  total  volume  at  these  operating  conditions. 

Rudinger's  analytical  results  (Eqns .  1.5,  1.6,  and  1.7)  for  the  isothermal 
case  can  be  used  to  calculate  nozzle  parameters.  Fig.  1.8  shows  typical 
results  for  gas  velocity,  dimensionless  nozzle  area,  and  the  particle  volume 
fraction  as  a  function  of  pressure  (and  hence  axial  distance)  through  the 
nozzle,  for  the  expansion  of  a  two-phase  mixture  of  air  (with  a  molecular 
weight  (MW)  of  28)  and  a  liquid  with  an  inlet  particle  volume  fraction  of  0.2 
and  an  inlet  temperature  of  260  C  (500  F).  As  seen  by  inspecting  Eq.  1.5, 
the  only  quantity  which  is  a  function  of  the  specific  liquid/gas  combination  is 
the  gas  constant.  Hence,  only  the  molecular  weight  of  the  gas  need  be  specified 
to  use  Eq.  1.5,  and  the  results  are  independent  of  the  liquid  phase.  If  steam 
were  the  gas  (MW  *  13)  at  the  same  pressure  level,  then  the  velocities  would  be 
higher,  and  the  required  flow  areas  through  the  nozzle  would  be  smaller  fractions 
of  the  inlet  area  to  pass  the  same  flow  rate,  as  deduced  from  Eq.  1.7.  In 
other  words,  a  nozzle  for  a  steam/liquid  flow  would  be  relatively  more  constricted 
at  the  throat  than  that  for  an  air/liquid  nozzle.  For  these  constant-fractional- 
lag  nozzles,  the  droplet  velocities  in  this  example  are  60  percent  of  the  local 
gas  velocities  in  the  nozzle. 

Equation  1.5  can  be  used  to  determine  the  effect  of  particle  volume  fraction  on 
the  gas  exit  velocity,  as  compared  to  the  case  when  e0  “  0  (the  all-gas  case), 
and  Fig.  1.9  shows  the  results.  For  values  of  eQ  less  than  10  percent, 
the  incxease  in  calculated  gas  velocity  is  less  than  4  percent  when  the  effect 
is  included.  At  larger  values  of  eQ,  however,  the  effect  is  not 
c.  igible.  For  example,  at  eQ  *  0.4,  the  increase  in  gas  velocity  is  12 
p;  -  nt  when  the  pressure  ratio  (PQ/P)  is  10. 

The  pressure  ratio  at  the  nozzle  throat  can  be  obtained  analytically  by 
determining  the  point  at  which  the  derivative  of  throat  area  with  respect  to 
pressure  ratio  goes  to  zero.  When  this  is  done,  the  effect  of  inlet  particle 
volume  fraction  (eQ)  on  the  throat  pressure  ratio  (Pt/PQ)  can  be  estimated 
with  the  isothermal  model,  and  Fig.  I. 10  shows  calculated  results  of  this 
effect.  For  the  all-gas  case  (eQ  ■  0),  the  throat  pressure  ratio  is  0.607, 
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and  as  the  inlet  particle  volume  fraction  is  increased,  the  throat  pressure 
ratio  decreases  from  the  all-gas  value.  For  models  which  neglect  the  particle 
volume  fraction,  the  predicted  throat  pressure  ratio  is  simply  0.607.  For 
example,  at  eQ  =  0.4,  the  throat  pressure  ratio  from  this  isothermal  model  is 
approximately  0.54. 

Calculated  results  from  the  homogeneous  and  the  separate-phase  models  can 
be  compared  with  each  other,  and  Fig.  I. 11  shows  calculated  results  for  gas  and 
liquid  velocities  as  a  function  of  exit  pressure  from  three  different  models 
for  a  two-phase  mixture  of  nitrogen  and  water  with  inlet  conditions  of  10.3  bar 
(150  psia)  and  10  C  (60  F)  expanding  through  a  converging-diverging  nozzle 
at  a  loading  ratio  of  40.  These  conditions  were  selected  because  Elliot  and 
Weinberg  (Ref.  8)  present  results  for  these  conditions  from  their  computerized 
comprehensive  model  which  includes  the  effects  of  velocity  slip,  temperature 
difference  between  the  phases,  and  droplet  breakup.  Results  for  the  mixture 
velocity  from  the  homogeneous  model  are  shown.  Because  the  value  of  the 
loading  ratio  is  40,  the  velocity  of  the  mixture  is  essentially  that  of  the 
liquid.  Also  shown  are  calculated  results  from  Rudinger's  isothermal  model,  at 
two  values  of  the  velocity  ratio  (V^/Vg):  0.71  and  1.0.  Rudinger's  model 
can  be  considered  as  partially  separated,  because  slip  is  included.  The  value 
of  0.71  was  selected  because  it  is  characteristic  of  the  velocity  ratio  over 
most  of  the  length  of  the  nozzle  of  Elliott  and  Weinberg.  The  calculated  gas 
and  liquid  velocities  from  Elliott  and  Weinberg's  model  are  shown  as  triangles. 
The  estimates  for  the  gas  and  liquid  velocities  from  Rudinger's  isothermal 
model  at  V^/Vg  *  0.71  coincide  closely  with  the  corresponding  results  from 
Elliott  and  Weinberg's  comprehensive  model,  and  the  prediction  from  the  homogen¬ 
eous  model  agrees  well  with  the  liquid  predictions  of  the  other  two  models. 

Results  calculated  from  the  homogeneous  model  and  Rudinger's  model  are 
shown  in  Fig.  1.12  for  expansion  of  a  steam/Dowtherm  A  mixture  with  typical 
inlet  conditions  and  a  loading  ratio  of  40.  (Predictions  from  Elliott  and 
Weinberg's  model  are  not  shown  because  a  listing  of  their  comprehensive  computer¬ 
ized  code  is  not  available  in  the  open  literature.)  For  these  conditions  at  an 
exit  pressure  of  P/PQ  =  0.2,  the  homogeneous  model  predicts  velocities 
about  16  percent  lower  than  the  liquid  velocity  from  the  isothermal  model  at  a 
slip  ratio  of  0.8. 

Figure  1.13  is  a  comparison  of  the  momentum  flux  at  the  nozzle  exit  as 
calculated  with  the  homogeneous  theory  (Eq.  1.4)  and  with  Rudinger's  isothermal 
theory.  The  total  momentum  of  the  two-phase  mixture  is  the  sum  of  the  momentums 
of  the  gas  and  liquid  streams,  and  the  momentum  flux  is  the  total  momentum 
divided  by  the  flow  area.  With  this  definition  and  the  equation  of  continuity, 
the  momentum  flux  can  be  shown  to  be: 
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Now,  ■  KVg,  and  the  gas  exit  velocity  for  the  isothermal  model  can  be 
calculated  with  Eq.  1.5  as  a  function  of  the  selected  operating  conditions.  As 
seen  in  Fig.  1.13,  the  iosthermal  model  predicts  higher  momentum  fluxes  at 
loading  ratios  larger  than  approximately  5  at  the  operating  conditions  typical 
for  the  marine  propulsion  application. 


An  estimate  of  the  effects  of  the  heat  transfer  rate  process  on  the 
temperature  difference  between  the  liquid  and  the  gas  as  the  gas  expands 
through  the  nozzle  can  be  made  using  a  simple  model.  Heat  transfer  is  a  rate 
process,  and  the  residence  time  determines  the  quantity  of  heat  transferred. 

For  heat  transfer  to  occur  between  the  liquid  and  the  gas,  there  must  be  a 
temperature  difference.  Even  if  the  liquid  and  the  gas  enter  the  nozzle  at  the 
same  temperature,  the  temperature  of  the  gas  will  drop  relative  to  the  liquid 
as  the  gas  expands,  thereby  creating  the  temperature  differential  which  is  the 
driving  force  for  heat  transfer.  To  make  numerical  estimates,  some  approxima¬ 
tions  can  be  made:  as  the  liquid  and  gas  flow  along  an  axial  increment  through 
the  nozzle,  assume  that  the  liquid  temperature  is  constant  over  this  increment, 
and  that  the  gas  temperature  drops  according  to  an  (adiabatic)  isentropic 
expansion.  Then,  for  this  small  axial  increment,  the  amount  of  heat  transfer 
from  the  gas  to  the  liquid  can  be  calculated,  and  the  temperatures  of  the  gas 
and  the  liquid  can  be  adjusted  at  the  end  of  this  increment  to  account  for  this 
heat  transfer.  With  these  assumptions  and  for  an  assumed  pressure  drop  over 
the  axial  increment,  the  gas  temperature  at  the  end  of  the  increment  can  be 
calculated  with  the  classical  isentropic  expression: 


Figure  1.14  depicts  the  model,  where  the  star  conditions  reflect  the 
temperatures  as  adjusted  at  the  end  of  the  increment  for  heat  transfer.  The 
heat  balance  on  a  liquid  droplet  i6: 


^drop  “  ■** 


(I. 11) 


where  m  is  the  mass  of  the  droplet,  C  is  its  specific  heat,  h  is  the  heat 
transfer  coefficient  between  the  droplet  and  the  gas,  A^  is  the  heat  transfer 
surface  area  of  the  droplet,  ATfc_g  is  the  average  temperature  differential 
between  the  liquid  and  the  gas,  and  &t  is  the  residence  time  of  the  droplet  in 
the  axial  increment  (At  *  aX/V^).  The  value  of  AT^_g  can  be  estimated  simply  as 


(1.12) 
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For  relative  droplet  Reynolds  numbers  in  the  range  of  25  <  ReRg^  <  100,000, 
the  heat  transfer  coefficient  can  be  estimated  with  the  relation  (as  recommended 
in  Ref.  19): 


M  -  0.37  Re0'6 

kg  *el 


(1.13) 


The  mass  of  the  droplet  is  P^'  Vol,  and  for  spherical  droplets,  Vol/A^  *  D.  Hence 

the  value  of  T*)out  can  be  calculated  from  Eq.  I. 11.  The  value  of  the 

gas  exit  temperature  Tg)  ,  adjusted  from  the  isentropic  value  to  account 

for  heat  transfer,  can  be  calculated  from  a  heat  balance  on  the  liquid  and  the 

gas: 

"£c  [Tt)ln  -  .  Vl>!  (Tg)o„t  -  T,)ootJ  (I-U) 

The  ratio  of  the  mass  of  liquid  droplets  (m^)  to  the  mass  of  gas  (mg)  in  the  axial 
increment  is  simply  equal  to  the  loading  ratio,  r.  Hence,  Eq.  1.14  reduces  to: 

Tg)out  *  Tg)out  +  r  |Ti)in  ~  T*)°utj  (1.15) 

and  the  adjusted  gas  exit  temperature  can  be  calculated.  This  simplified  model 
was  used  to  estimate  the  effect  of  the  heat  transfer  rate  process  on  the 
temperature  difference  between  the  liquid  and  gas  at  typical  operating  conditions. 
Fig.  1.15  shows  such  estimates  for  mixtures  of  steam  and  Dowtherm  A  at 
an  inlet  temperature  of  260  C  (500  F)  as  a  function  of  loading  ratio,  droplet 
diameter,  and  velocity  level  at  a  slip  ratio  (V^/Vg)  of  0.8.  The  temperature 
of  247.1  C  (476.9  F)  occurs  when  steam  is  expanded  over  a  pressure  of  P/PQ  =  0.9. 
Three  factors  can  cause  the  gas  temperature  in  the  nozzle  to  approach  that  of 
the  liquid:  first,  a  decrease  in  drop  diameter,  because  the  heat  transfer 
coefficient  and  the  volumetric  heat  transfer  area  both  increase;  second,  an 
increase  in  loading  ratio,  because  the  amount  of  gas  relative  to  the  amount  of 
liquid  is  reduced;  and  third,  a  decrease  in  velocity,  because  the  residence 
time  increases  faster  than  the  heat  transfer  coefficient  decreases.  Actually, 
the  model  is  over-simplified  in  that  it  does  not  predict  the  asymptote  of 
960  R  as  the  loading  ratio  is  increased  and  the  asymptote  must  be  faired 
in,  but  the  model  does  indicate  the  conditions  at  which  the  heat  transfer  rate 
can  have  an  effect  on  nozzle  performance.  Fig.  1.16  shows  calculations  from 
the  model  when  the  inlet  temperature  of  the  Dowtherm  A  is  56  C  (100  F)  above 
the  steam  inlet  temperature,  and  the  same  trends  occur  as  in  Fig.  1.15.  this 
simplified  heat  transfer  model  indicates  the  conditions  under  which  the  flow  is 
likely  to  be  isothermal  and  hence  under  which  the  isothermal  theory  is  valid. 
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Gas  dynamics  aspects  of  nozzle  flow  were  examined  in  terms  of  the  speed  of 
sound  in  a  two-phase  mixture  and  choked  flow  at  the  nozzle  throat.  The  change 
in  operating  conditions  at  one  point  in  a  flow  system  is  telegraphed  to  other 
points  by  pressure  waves  which  travel  at  the  speed  of  sound  and  result  in  the 
necessary  flow  adjustments.  If  thermodynamic  equilibrium  (temperature,  pressure, 
and  kinetic  equilibrium)  between  phases  is  assumed  for  a  homogeneous  two-phase 
mixture,  then  the  sonic  velocity  can  be  defined  as  the  square  root  of  the 
isentropic  gradient  of  the  pressure  with  respect  to  density,  as  in  the  case  of 
the  classical  gas  dynamics  theory.  With  this  definition,  the  two-phase  choking 
velocity  is  always  lower  than  the  choking  velocity  of  the  pure  saturated  vapor 
at  the  same  temperature  and  pressure.  Rudinger  (Ref.  20)  has  shown  that  the 
equilibrium  speed  of  sound  a£  in  the  mixture  relative  to  that  of  the  pure 
vapor  aQ  is : 


where  the  last  term  accounts  for  the  finite  droplet  volume  fraction  at  high 
values  of  the  loading  ratio.  Fig.  1.17  shows  typical  results  from  Eq. 

1.16  for  a  mixture  of  saturated  steam  and  Dowtherm  A  at  260  c  (500  F). 

The  two-phase-mixture  speed  of  sound  reaches  a  minimum  of  0.29  aQ,  at  a  loading 
ratio  of  approximately  20.  The  value  of  aQ  at  260  C  (500  F)  is  566 
m/sec  (1857  ft/sec),  so  the  choked  velocity  of  a  homogeneous  mixture  is  164 
m/sec  (539  ft/sec)  in  the  nozzle  throat.  As  seen  earlier  in  Fig.  1.7,  the 
inlet  particle  volume  fraction  for  this  mixture  is  expected  to  be  between  0.4 
and  0.48  at  a  loading  ratio  of  20.  Hence,  from  Fig.  1.10,  the  throat  pressure 
ratio  (Pt/PQ)  is  between  0.52  and  0.54  for  this  range  of  inlet  particle 
volume  fractions.  From  Fig.  1.17,  for  a  pressure  ratio  in  this  range,  the 
isothermal  model  predicts  a  gas  velocity  of  between  137  and  140  m/sec  (450 
and  460  ft/sec),  at  a  slip  ratio  (V£/Vg)  of  0.8.  At  this  same  pressure 
ratio,  the  homogeneous  model  predicts  a  velocity  of  between  85  and  88  m/sec 
(280  and  290  ft/sec).  These  estimates  suggest  that  the  mixture  velocity  at 
the  nozzle  throat  at  these  operating  conditions  is  below  the  choked  flow 
velocity. 

Comfort  and  Crowe  (Ref.  21)  investigated  the  effect  of  droplet  size  on 
shock  characteristics  of  a  steam/water  flow  in  a  converging-diverging  nozzle. 

As  the  droplet  size  is  reduced,  the  mixture  behaves  as  a  continuum  and  sharp 
velocity  discontinuities  occur  at  velocities  above  the  two-phase  choking  velocity 
but  below  the  vapor  sonic  velocity.  Analytically,  for  their  conditions  and 
droplet  diameters  below  about  2  microns  (7.8  x  10“  ^  in.),  then  the  mixture 


R80-954624-2 


behaves  as  a  continuum,  and  the  shock  conditions  can  be  carried  upstream 
because  the  sonic  velocity  is  higher  than  the  actual  gas  velocity.  For 
larger  droplet  diameters,  the  change  in  velocity  is  very  gradual,  without 
discontinuity  as  the  mixture  no  longer  behaves  as  a  continuum.  Because  the 
mass  fraction  of  the  dispersed  liquid  droplets  strongly  affects  the  choking 
velocity,  relatively  low-velocity  flows  can  be  supersonic  with  respect  to  the 
homogeneous  two-phase  choking  velocity.  At  another  operating  condition, 
where  the  gas  velocity  is  above  the  vapor  sonic  velocity.  Comfort  and  Crowe 
calculated  that  the  maximum  isentropic  velocity  (of  the  vapor)  in  the  diverging 
section  of  their  nozzle  with  steam/water  flow  would  be  approximately  Mach  1.2 
relative  to  the  vapor  sonic  velocity  and  Mach  3.3  relative  to  the  two-phase 
choking  velocity  for  a  homogeneous-mixture  model.  Here  the  droplet  size  has 
little  effect  on  the  sharpness  of  the  shock,  because  the  shock  is  not  felt 
upstream  when  the  gas  velocity  is  above  the  vapor  sonic  velocity.  Hence, 
unless  the  droplets  are  very  small  or  unless  the  gas  velocity  is  above  the 
vapor  sonic  velocity,  the  flow  disturbance  is  very  gradual  and  the  sharp 
discontinuities  of-  classical  gas  dynamics  are  not  expected.  Rudinger  (Ref.  22) 
has  shown  that  the  distance  that  a  droplet  travels  while  achieving  velocity  and 
temperature  equilibrium  with  the  gas  is  orders  of  magnitude  larger  than  the 
distance  for  shock  transition  in  a  pure  gas,  and  this  is  the  reason  why  flow 
disturbances  attributable  to  shock  conditions  in  a  two-phase  mixture  can  often 
be  just  gradual  changes  in  flow  conditions. 

1.3.4  Correlation  of  Nozzle  Flow  Domains 


The  various  nozzle  flow  domains  were  correlated  and  their  physical 
characteristics  were  identified  in  terms  of  the  governing  flow  parameters. 

For  many  of  the  two-phase-nozzle  flow  conditions,  the  flow  regime  will  be 
that  of  dispersed  flow,  where  the  liquid  phase  is  in  the  form  of  uniformly 
dispersed,  small-diameter  droplets.  Dispersed  flow  through  the  nozzle  is 
desirable,  because  uniformly  dispersed  small  droplets  result  in  enhanced  heat 
transfer  between  the  liquid  and  the  gas  phases  as  the  droplet  heat  transfer 
coefficient  and  the  volumetric  heat  transfer  surface  area  both  increase  as 
the  droplet  diameter  is  reduced.  The  assumption  of  a  homogeneous  flow  model 
implies  small  droplets  uniformly  dispersed  over  the  croBS-sect ion  of  the  flow 
area  and  in  thermal  equilibrium  with  the  gas.  In  the  limit,  if  the  diameter 
of  the  droplets  and  their  average  spacing  are  both  small  as  compared  to  the 
nozzle  throat  diameter,  then  the  mixture  begins  to  behave  as  a  continuum.  At 
a  given  loading,  as  the  number  of  droplets  is  reduced,  then  the  flow  field 
becomes  less  uniform,  and  the  limit  is  one  droplet  where  it  is  apparent  that 
the  homogeneous  theory  no  longer  applies.  The  assumption  of  homogeneous  flow 
conditions  is  sometimes  made  even  though  the  flow  is  not  homogeneous,  because 
useful  results  can  be  obtained  from  the  homogeneous  flow  model.  The  approxi¬ 
mation  of  the  homogeneous  model  can  only  be  applied  to  flow  situations  where 
the  two  phases  are  well  mixed,  and  this  condition  excludes  such  flow  regimes 
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as  slug,  annular,  or  stratified  flows.  And  the  presence  of  any  one  of  these 
flow  regimes  would  suggest  poor  nozzle  design  in  the  first  place,  as  their 
presence  would  negate  the  advantages  of  dispersed  flow.  Based  on  the  work  by 
Comfort  and  Crowe  (Ref.  21)  that  was  discussed  above  in  Section  1.3.3, 
however,  the  droplet  size  at  which  the  mixture  really  behaves  as  a  continuum 
in  terms  of  shock  characteristics  is  very  small,  on  the  order  of  less  than  3 
microns  (2  x  10-^  in.)  in  diameter.  With  reference  to  the  Weber  number 
droplet  breakup  criterion  shown  in  Fig.  1.4,  this  size  droplet  would  be 
expected  at  gas  velocities  higher  than  approximately  180  m/sec  (600  ft/sec), 
and  this  velocity  would  probably  exist  near  the  exit  of  many  two-phase 
nozzles.  Even  though  truly  homogeneous  conditions  may  not  obtain  throughout 
the  nozzle,  the  temperature  and/or  velocity  of  the  two  phases  may  be  in 
sufficient  equilibrium  so  that  results  from  the  homogeneous  model  could  still 
be  useful  for  examining  many  nozzle  applications. 

As  indicated  by  the  results  obtained  from  the  simplified  heat  transfer 
model  of  Fig.  1.14,  the  two-phase  flow  becomes  increasingly  isothermal  as  the 
droplet  size  is  reduced  (because  the  heat  transfer  rate  is  increased)  and  as 
the  loading  ratio  is  increased  (because  the  increased  amount  of  liquid 
relative  to  the  gas  is  a  source  of  heat).  As  typified  by  the  results  in  Fig. 
1.2,  the  homogeneous  theory,  with  infinite  heat  transfer,  indicates  the 
values  of  the  loading  ratio  at  which  isothermal  conditions  can  be  approached 
(depending  on  the  droplet  diameter),  but  the  homogeneous  model  does  not 
account  for  the  fact  that  heat  transfer  is  a  rate  process  which  requires 
residence  time.  For  flow  conditions  which  do  lead  to  isothermal  flow, 
however,  Rudinger's  isothermal  model  is  applicable  and  still  allows  for  slip 
between  the  liquid  and  the  gas. 


An  indication  of  the  degree  of  velocity  slip  can  be  obtained  by  looking 
at  the  relaxation  time  for  the  particle  velocity,  Ty,  which  is  a  qualitative 
measure  of  the  time  to  achieve  velocity  equilibrium  between  the  gas  and  the 
droplet.  Consider  the  motion  of  a  single  droplet  in  a  gas  flow,  and  asstime 
that  Stokes  drag  is  valid.  Equate  the  drag  of  the  gas  with  constant  velocity, 
Vg,  to  the  droplet  acceleration: 


^  (Vg  -  V«) 


(1.17) 


The  solution  to  this  simplified  equation  is: 


VV*)c 


e“t/T 


(1.18) 
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where  Vl>0  is  the  initial  velocity  at  time  zero,  and  ty  is  the  relaxation  time 
for  the  particle  velocity: 


t 


v 


(1.19) 


Now,  the  residence  time  of  the  mixture  with  average  gas  velocity,  Vg)ave»  in 
the  nozzle  of  length  L  is:  treg  ■  L/Vg^ave 

Hence  a  dimensionless  residence  time  can  be  formed  by  the  ratio: 


'ave 


(1.20) 


When  the  residence  time  is  large  compared  to  the  velocity  relaxation  time,  then 
velocity  equilibrium  between  the  gas  the  the  liquid  can  be  expected.  Fig. 

1.18  shows  Eq.  1.20  plotted  for  typical  conditions  of  a  steam/Dowtherm  A 
mixture.  Several  trends  can  be  observed:  the  direction  of  increasing 
velocity  equilibrium  is  smaller  droplet  diameter,  lower  gas  velocities,  and 
longer  nozzles,  with  droplet  diameter  being  the  single  most  important  factor. 
For  example,  with  a  0.15-m  (6-inch)  long  nozzle  and  an  average  gas  velocity 
of  152  m/sec  (500  ft/sec),  the  droplet  diameter  would  have  to  be  less  than 
12.5  microns  (0.0005  inches)  to  approach  conditions  of  velocity  equilibrium 
between  the  gas  and  the  droplet. 


1.4  Formulation  of  Working  Theory  For  Two-Phase  Nozzles 

The  objective  of  this  task  was  to  formulate  a  working  theory  of  two-phase 
nozzle  flow  based  on  considerations  of  the  governing  flow  parameters  and  the 
flow  domains.  The  theory  had  to  be  in  a  form  suitable  for  parametric  analyses 
and  optimization  of  two-phase  flow  nozzles.  Two  simplified  models  have  been 
developed  and  computerized:  an  Isothermal  Model  which  is  valid  only  for  certain 
operating  conditions;  and  a  Model  with  Droplet  Heat  Transfer  which  is  a  more 
general  model  and  is  valid  for  a  wide  range  of  operating  conditions  since  it 
includes  the  effects  of  heat  transfer  between  the  liquid  droplets  and  the  gas 
phase.  Both  models  include  the  interphase  effects  of  velocity  slip  and  droplet 
breakup. 

Both  models  assume  a  two-component,  two-phase  mixture,  where  the  liquid 
phase  is  in  the  form  of  droplets.  One-dimensional  flow,  with  variables 
changing  only  in  the  axial  flow  direction  of  the  nozzle,  was  selected  so 
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both  the  model  and  the  results  will  bz  tractable  in  terms  of  parametric 
analysis  and  optimization  of  nozzles.  A  two-component  mixture  was  chosen 
for  the  marine  propulsion  application  because  the  expansion  of  a  two-component 
system  can  be  nearly  isothermal  (depending  on  the  loading  ratio)  when  compared 
with  a  one-component  system  which  expands  along  its  saturation  vapor  line.  As 
discussed  in  Section  1.2.1,  more  area  is  enclosed  by  the  integral  of  TdS  (Fig. 

1.3)  when  the  expansion  is  nearly  isothermal  than  when  an  adiabatic,  one-component 
expansion  is  experienced. 

Numerical  marching-type  solutions  have  been  selected  for  both  models 
because  closed-form  analytical  solutions  are  possible  only  for  certain  simplified 
situations  as  discussed  in  Section  1.3.  In  computing  the  solutions,  the  nozzle 
is  divided  into  numerous  small  axial  segments  of  equal  length,  and  the  flow 
conditions  at  the  exit  of  each  segment  are  calculated  from  the  governing  flow 
equations  and  the  upstream  flow  conditions  at  the  inlet  of  that  axial  segment. 

The  exit  flow  conditions  of  a  given  segment  then  become  the  inlet  flow  conditions 
for  the  adjacent  downstream  segment,  hence  the  description,  "marching-type" 
solution,  as  the  calculations  proceed  from  one  segment  to  the  next  downstream 
segment.  A  six-inch-long  nozzle,  for  example,  could  be  divided  into  1000 
segments,  each  0.006  in.  long  (which  could  be  larger  than  the  droplet  diameter 
for  many  flow  conditions).  Care  must  be  shown  to  select  enough  segments  so 
that  the  equations  are  numerically  stable. 

In  the  Model  with  Droplet  Heat  Transfer,  the  variations  of  the  thermo- 
physical  properties  of  the  gas  and  the  liquid  with  temperature  and/or  pressure 
along  the  nozzle  are  included  in  the  marching-type  solution  by  recalculating 
the  value  of  the  property  for  the  local  temperature  or  pressure  at  each  segment. 

In  the  Isothermal  Model,  the  thermophysical  properties  are  not  temperature- 
dependent  because  the  flow  is  assumed  as  isothermal. 

1.4.1  Isothermal  Model 

In  the  Isothermal  Model,  the  assumption  is  made  that  the  liquid  and  gas 
temperatures  are  equal  and  constant  throughout  the  nozzle  (T^  ■  Tg  *  T^n).  This 
assumption  implies  that  the  heat  transfer  rate  between  the  liquid  droplets  and 
the  gas  phase  is  very  high  and  that  the  liquid  droplets  act  as  a  heat  source  so 
that  the  gas  temperature  does  not  decrease  as  the  gas  expands  through  the  nozzle. 
This  assumption  means  that  the  Isothermal  Model  is  suitable  only  for  certain 
operating  conditions,  such  as  very  small  droplet  size  and  high  values  of  the 
loading  ratio  (m^/mg).  For  the  isothermal  situation  to  occur  even 
in  theory,  it  would  require  either  an  internal  heat  source,  or  an  external  heat 
input,  for  the  expansion  of  the  two-phase  mixture  in  the  nozzle  to  be  isothermal. 


kv> 
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In  addition  to  the  Isothermal  Model  being  valid  for  certain  operating 
conditions,  this  model  serves  as  a  useful  adjunct  to  a  model  that  does  account 
for  finite  droplet  heat  transfer,  that  is,  it  is  useful  because  the  calculated 
results  from  its  momentum  equations  can  be  compared  with  the  results  from  the 
analytical  isothermal  model  of  Rudinger  (Ref.  4). 

The  following  four  basic  equations  describe  the  isothermal  one-dimensional, 
two-phase  nozzle  flow: 


Cont inuity 
Momentum 
Droplet  Drag 
Droplet  Breakup 


These  equations  describe  the  flow  conditions  in  and  across  each  segment  of 
finite  length,  and  each  has  been  transformed  into  a  form  suitable  for 
computerized  numerical  solutions. 


The  Isothermal  Model  makes  allowance  for  velocity  slip  between  the 
liquid  droplets  and  the  gas  and  also  includes  the  effects  of  droplet  drag 
and  droplet  breakup.  The  nozzle  inlet  conditions  are  specified  as  input 
to  the  model,  and  the  pressure  distribution  along  the  axial  length  of  the 
nozzle  is  also  required. 


1 . 4 . 1 . 1  Coti t^nui  t  i°IL 


At  a  given  nozzle  cross-section,  the  total  flow  area  (A)  can  be 
considered  as  the  sum  of  the  liquid  flow  area  (A£)  and  the  gas  flow 
area  (A^): 


A  ■  A„  +  A 
1  g 


(1.21) 


For  one-dimensional  flow  through  a  cross-section,  the  continuity  equations 
for  each  of  the  two  components  can  be  written  as: 


m 


p  V  A 
g  g  g 


p  V  A 
l  i  l 


(1.22) 


These  forms  of  the  continuity  equation  imply  the  assumptions  that  the  gas 
does  not  dissolve  in  the  liquid  and  that  the  vapor  pressure  of  the  liquid 
is  so  low  that  the  liquid  does  not  vaporize. 
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Recalling  Che  definition  of  Che  loading  ratio  (r  s  m^/nig), 
then  the  expressions  in  Eq.  1.22  can  be  substituted  into  Eq.  1.21  to  produce 


A  ■  m 

g 


(1.23) 


This  is  the  form  of  the  continuity  equation  that  is  used  with  the  momentum 
equation. 


1 . 4 . 1 . 2  Home n t um_E^uat ion 

A  differential  force  balance  on  a  control  volume  with  cross-section  A 
and  infinitesimal  axial  length  dx  describes  the  conversion  of  pressure 
drop  (dP)  into  momentum: 

dM  *  -A  dP  (1.24) 

where  the  total  momentum  (m)  of  the  two-phase  mixture  is  equal  to 
the  sum  of  the  momentum  of  each  of  the  two  streams: 


ft 


mV  +  m  V 

g  g  t.  e. 


The  mass-weighted  mean  velocity  (V)  is  defined  as: 


V  H 


V  +  rV 

_8 _ L 

1  +  r 


(1.25) 


(1.26) 


Hence,  in  terms  of  the  weighted  velocity,  Eq.  1.26  for  the  total  momentum 
becomes : 


ft  - 


V 


(1.27) 


Taking  the  differential  of  Eq.  1.27,  and  recalling  that  (mg  +  m&) 
is  constant  by  conservation  of  mass,  results  in: 


dft  - 


(1.28) 
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Then,  substituting  Eqns.  1.28  and  1.23  into  the  basic  force  balance  of 
Eq.  1.24  yields  a  differential  expression  for  the  weighted  velocity: 

Actually,  V  is  a  function  of  the  phase  velocities,  V  and  V^,  and 
expressions  which  relate  these  three  velocities  are  developed  below. 

The  local  velocity  ratio  (K)  is  defined  as:  K  =  V^/V  .  Combining 
this  definition  with  the  mean-velocity  definition  (Eq.  1.26)  produces 
expressions  for  the  local  phase  velocities,  Vt  and  Vg,  in  terms  of  the 
local  velocity  ratio  and  the  local  mean  velocity: 


(1.30) 


Using  these  two  expressions,  Eq.  1.29  can  be  rewritten  as: 


(1.31) 


(1.32) 


Equation  1.32  can  now  be  transformed  into  a  form  for  numerical  integration: 

P+AP 

«•*> -*’*-/  (rh 

p 

where  the  perfect  gas  law  is  assumed  for  the  density  of  the  gas  phase 
(Pg  *  PMg/RTg),  and  where  AP  is  the  pressure  change  over  the  length 
of  the  segment  AX  as  the  two-phase  mixture  flows  from  position  X  to 
downstream  position,  X  ♦  AX.  Mg  is  the  molecular  weight  of  the  gas 
and  R  is  the  universal  gas  constant.  If  the  assumption  is  made  that 
the  quantities  Tg,  G,  L,  and  p(  vary  only  slowly  over  the  small 
pressure  increment  AP,  then  the  values  of  these  quantities  can  be 
assumed  to  be  constant  over  the  small  segment  and  can  be  assigned  the 
numerical  value  that  exists  at  the  inlet  of  each  segment.  The  values 
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of  these  quantities  are  then  recalculated  for  each  segment.  With  this 
assumption,  Eq.  1.33  can  be  integrated  to: 


This  is  the  form  of  the  momentum  equation  that  is  used  in  the  computerized 
numerical  solution. 


1 . 4 . 1 . 3  Dropl_e t_  J)r £g_Equat icm 

Droplet  drag  provides  the  mechanism  by  which  the  gas-phase  momentum, 
which  increases  during  expansion,  is  transferred  to  the  liquid  droplets. 
Accordingly,  this  momentum  transfer  produces  droplet  acceleration.  Using 
a  force  balance  on  the  droplet,  as  shown  schematically  in  Fig.  1.19,  the 
droplet  acceleration  force  is  equal  to  the  drag  force  less  the  force  caused  by 
the  pressure  gradient  of  the  gas  phase: 


F  -  F  =  F 

Drag  Press  Accel 


(1.35) 


Because  the  pressure  gradient  is  negative,  the  pressure  is  lower  on  the 
downstream  side  of  the  droplet  than  on  the  upstream  side,  and  the  overall 
effect  of  the  pressure-gradient  force  is,  therefore,  to  accelerate  the  droplet. 

The  drag  force  on  the  droplet  is  represented  by  the  expression: 


F  «  —  p  V  !v  Ic 
Drag  2  Pg  Slip 1 VSlip 1  S>  4 


(1.36) 


where  Vg££p  is  the  velocity  slip  between  the  two  phases  (Vg£-  ■  Vg  -  V^)  and 

the  absolute  symbol  around  Vgj-  accounts  for  the  direction  of  the  drag  force. 

If  Vt  is  larger  than  Vg,  then  tne  drag  force  on  the  droplet  would  be 
negative  (and  the  direction  of  the  force  is  properly  predicted  by  Eq.  1.36).  A 
convenient  definition  (Ref.  8)  is  the  velocity  slip  relative  to  the  mass-weighted 
mean  velocity: 


S 


(1.37) 
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The  pressure-gradient  force  on  a  small  droplet  in  a  pressure  gradient 
is  represented  by  the  expression  (Ref.  23): 


F 


Press 


dP 
'D  dX 


V 


(1.38) 


2 

where  is  the  volume  of  the  droplet.  For  spherical  droplets*  Vp  =  xD  /6, 
and  Eq.  1.37  becomes: 


F 


Press 


ttD3  dP 
6  dX 


(I. 38a) 


Even  if  the  velocity  slip  (Vg^£p)  is  zero  and  the  drag  force  is  therefore 
zero,  the  pressure  gradient  caused  by  the  gas  phase  still  exists  and  exerts 
a  force  on  the  droplet. 


The  acceleration  force  on  the  droplet  is: 


Accel 


ma  = 


(1.39) 


Combining  the  three  force  equations  into  the  force  balance  of  Eq.  1.35 
and  using  the  convenient  definition  of  Eq.  1.37,  the  following  equation 
results : 


% 


3p  Isis  V2  C  dX 
g  P 

4  ° 


dP 


P  V 


(1.40) 


The  definitions  of  the  mean  velocity  (Eq.  1.26)  and  the  slip  (Eq.  1.37)  can 
be  combined  to  eliminate  the  explicit  variable  V  ,  thereby  obtaining  an 
expression  for  in  terms  of  V  and  s: 


V 


l 


(1.41) 


This  expression  is  consistent  with  Eq.  1.30  for  in  terms  of  V  and 

K.  For  the  cases  considered  here,  the  assumption  is  made  that  no  mass  transfer 

occurs  between  the  phases.  Hence  the  loading  ratio  (r)  is  constant  throughout 
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the  nozzle.  Noting  this  and  taking  the  first  differential  of  Eq.  1.41  results 
in: 


dV* 


V 

(1  +  r) 


dS 


(1.42) 


Then,  Eqns .  1.42  and  1.40  can  be  combined  to  eliminate  the  term,  dV^: 


dS 


ho  |s|s  V2  cD  dx 

\  §4pAD 


(1.43) 


Recalling  that  =  LV  from  Eq.  1.30  and  noting  that  dV^  =  2VdV 

(as  ^  is  the  way  that  V  appears  in  Eq.  1.34  for  the  momentum) , then  Eq.  1.43 

takes  the  following  form  for  an  axial  segment  of  finite  length  AX: 


AS  = 


1  +  r 
2  V2 


1  +  r, 


3p  |S  S  C  (1  +  r)  AX 

-g  _ 5 _ 

4  P gLD 


O-  +  r)  AP 

+  _2 

p^V 


(1.44) 


This  equation  is  valid  for  numerical  analysis  if  the  axial  segments  are 
short  enough  so  that  the  equation  can  be  linearized  over  the  segment,  and 
quantities  can  be  evaluated  at  their  inlet  values  (except  for  differential 
quantities).  This  equation  represents  the  differential  change  of  slip  that 
occurs  as  momentum  is  transferred  from  the  gas  to  the  liquid  by  means  of  the 
interphase  coupling  from  the  drag  force.  The  drag  force  and  the  buoyancy  force 
caused  by  the  pressure  gradient  across  the  droplet  are  balanced  by  the  droplet 
acceleration  force. 


1.4. 1.4  Droplet  Breakup 


The  fourth  relationship  required  to  describe  the  two-phase  flow  in  the 
Isothermal  Model  is  the  droplet  breakup  criterion.  Droplet  breakup  in  a 
two-phase  mixture  is  governed  primarily  by  the  ratio  of  two  forces:  the  aero¬ 
dynamic  pressure  forces  which  are  trying  to  deform  the  droplet,  and  the  surface 
tension  forces  which  are  trying  to  maintain  the  droplet  in  spherical  form.  The 
ratio  of  these  two  forces  is  known  as  the  Weber  Number,  We,  defined  as: 


We 


2 

VsliP 

2o 


D 


(1.45) 
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At  high  values  of  the  Weber  Number,  the  aerodynamic  forces  are  large,  and 

the  droplets  will  deform  and  fragment  into  smaller  droplets.  At  low  values  of 

the  Weber  Number,  the  surface  tension  forces  dominate  and  the  drops  remain 

nearly  spherical.  The  critical  value  of  the  Weber  Number,  We)(jrj.t»  i® 

experimentally  determined  and  represents  the  maximum  droplet  size  in  terms  of 

droplet  diameter,  D  ,  that  can  exist  at  the  local  flow  conditions  in  a 
’  max’ 

two-phase  flow: 


D 

max 


V 


2 

slip 


(1.46) 


The  generally  accepted  value  of  We)^r£t  is  6.0  (Ref.  23),  within  a  factor  of 
about  2  (for  example,  see  Ref.  24).  However,  recent  droplet  size  distribution 
measurements  by  Alger  (Kef.  13)  with  steam-water  nozzles  suggest  that  the 
appropriate  value  of  We)gr£t  should  be  approximately  1.0.  Regions  of  high 
acceleration  in  the  nozzle  could  produce  greater  droplet  deformation  and  hence 
earlier  breakup  as  compared  to  nearly  spherical  droplets  in  less  accelerated 
regions.  For  the  smaller  droplets  predicted  which  correspond  to  a  smaller 
value  of  We)crit»  the  interphase  coupling  of  momentum  (and  also  heat 
transfer  in  the  Model  with  Droplet  Heat  Transfer,  to  be  discussed  in 
Section  1.42)  will  be  enhanced. 


1 . 4 . 1 . 5  ija  n_T  ec  hn_i  q  uie 

The  four  basic  equations  used  in  the  Isothermal  Model  to  solve  the 
two-phase  flow  presented  earlier  include:  continuity  (Eq.  1.23);  momentum 
(Eq.  1.34);  droplet  drag  (Eq.  1.44);  and  droplet  breakup  (Eq.  1.46).  The 
nozzle  inlet  flow  conditions  (gas  and  liquid  velocities,  loading  ratio, 
temperature  and  pressure  of  the  mixture  and  the  pressure  drop  along  the 
nozzle  must  be  specified,  as  well  as  should  the  thermophysical  properties 
of  the  gas  (molecular  weight  and  viscosity)  and  of  the  liquid  (density  and 
surface  tension). 


The  solution  technique  is  a  marching-type  solution,  whereby  the 
governing  flow  equations  are  solved  for  the  first  segment,  and  the 
calculated  results  from  the  first  segment  are  then  used  as  input  for 
calculating  the  solution  for  the  second  segment,  and  so  on  through  the 
nozzle.  Basically,  the  quantity  is  calculated  for  the  segment 
using  Eq.  1.34  and  the  specified  pressure  distribution.  Hence,  the  value 
of  V  at  the  exit  of  the  segment  can  be  calculated: 


(1.47) 
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The  change  in  the  slip  ratio,  AS,  over  the  segment  can  then  be  calculated 
using  Eq.  1.44,  and  the  exit  value  of  the  slip  ratio  can  be  calculated  from: 


S) 


Exit 


S) 


Inlet 


+  AS 


(1.48) 


Using  the  exit  values  of  the  slip  ratio  and  the  weighted  velocity,  the  value 
of  the  liquid  velocity,  Vt,  at  the  exit  of  the  segment  can  be  calculated 
using  Eq.  1.41.  The  exit  value  of  the  gas  velocity  can  then  be  calculated  from 
the  definition  of  the  slip  ratio  (Eq.t1.37),  and  the  nozzle  area  can  be  calcu¬ 
lated  as  a  function  of  axial  distance  .using  Eq.  1.23  once  V  and  Vo  are 
known.  g 

Results  which  were  calculated  using  the  Isothermal  Model  for  selected 
operating  conditions  are  presented  and  discussed  later  in  this  section. 

1.4.2  Model  With  Droplet  Heat  Transfer 

The  main  difference  between  the  Isothermal  Model  discussed  above  and 
the  Model  with  Droplet  Heat  Transfer  is  that  the  latter  includes  the  effect  of 
heat  transfer  by  allowing  a  varying  temperature  difference  between  the  liquid 
droplets  and  the  gas  and  an  axial  temperature  decrease  in  the  liquid  and  the 
gas.  In  addition  to  the  four  basic  equations  cited  for  the  Isothermal  Model, 
two  additional  governing  equations,  the  energy  equation  and  an  equation  for 
heat  transfer  from  the  liquid  to  the  gas,  are  required  so  the  liquid  and  gas 
temperatures  of  the  two-phase  flow  can  be  calculated.  The  basic  equations 
required  are : 


Continuity 
Momentum 
Droplet  Drag 
Droplet  Breakup 

Energy 

Heat  Transfer 

The  first  four  equations  are  used  in  the  same  form  as  are  presented  for  the 
Isothermal  Model.  In  the  model  with  heat  transfer,  however,  the  thermo¬ 
physical  properties  are  temperature-dependent,  including  the  temperature- 
dependence  of  the  gas  density  based  on  the  ideal  gas  equation  that  was  used 
in  the  formulation  of  the  momentum  equation  (Section  I. 4. 1.2). 
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1.4. 2.1  Energy  Equation 


The  energy  equation  of  the  mixture  states  that  the  net  energy  flux 
into  a  volume  element  is  zero,  because  the  two-phase  flow  in  the  nozzle  is 
steady  state  and  assumed  to  be  adiabatic.  Hence,  the  total  enthalpy  of  the 
mixture  is  constant: 


V 

g 


dT 


+  C  — 1 B  +  r 

Pg  <*X 


V 


dv„  dT„ 

— -  +  C  — - 
i  dX  dX 


+ 


1_  dP 
p£  dX 


(1.49) 


The  first  two  terms  represent  the  gas  enthalpy,  and  the  other  three  terms 
represent  the  liquid  enthalpy.  Equation  1.49  applies  throughout  the  nozzle, 
and  the  integrated  form  of  the  energy  equation  is: 
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where  the  Constant  term  is  determined  from  nozzle  inlet  conditions.  The 
energy  equation  is  an  overall  balance  on  the  energy  of  the  two-phase  mixture. 
However,  the  temperatures  of  the  gas  and  of  the  liquid  can  differ,  yet  satisfying 
the  above  energy  equation. 


1 . 4 . 2 . 2  He  a  t__Tr ans  fejr  Equation 

As  the  gas  phase  of  the  two-phase  mixture  expands  through  the  nozzle,  the 
temperature  of  the  gas  will  decrease.  It  is  assumed  that  the  temperature  of 
the  liquid  phase  is  always  equal  to  or  greater  than  that  of  the  gas  phase  at 
the  nozzle  inlet.  In  this  case,  as  the  gas  temperature  drops,  the  temperature 
differential  between  the  liquid  and  the  gas  will  increase,  and  heat  will  be 
transferred  from  the  liquid  to  the  gas  at  an  increased  rate,  other  conditions 
being  equal.  A  heat  balance  on  the  liquid  droplet  describes  this  phenomenon: 
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(1.51) 


In  this  case,  the  rate  of  change  of  the  enthalpy  of  the  droplet  is  equal 
to  the  rate  of  heat  transfer  from  the  droplet  to  the  surrounding  gas.  The 
heat  transfer  coefficient,  h,  is  obtained  from  experimental  data.  For 
small  segments,  the  time  differential  can  be  approximated  as:  bt  *  &X/Vt ,  and 
for  spherical  droplets,  the  volumetric  density  of  heat  transfer  area  is: 
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Hence,  Eq.  1.51  can  be  solved  for  the  temperature  change  of  a  liquid  droplet 
in  the  axial  segment  of  length  AX: 
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Equation  1.53  is  then  used  to  account  for  droplet  heat  transfer  effects. 

I. 4. 2. 3  Solution  Technique 

The  basic  solution  technique  for  the  Model  with  Droplet  Heat  Transfer 
is  very  similar  to  that  discussed  in  Section  I. 4. 1.5  for  the  Isothermal 
Model,  except  that  the  energy  equation  (Eq.  1.50)  and  the  heat  transfer  equation 
(Eq.  1.53)  must  now  be  included. 

Again,  a  marching-type  solution  is  used.  Based  on  conditions  at  the  inlet 
to  the  segment,  the  temperature  change  of  the  liquid  droplets  is  calculated  by 
using  Eq.  1.53.  The  heat  transfer  coefficient  is  calculated  with  correlations 
for  droplet  heat  transfer  based  on  the  Reynold's  Number.  The  liquid  temperature 
at  the  exit  of  the  segment  is  then  known  from  the  calculated  change  in  tempera¬ 
ture  and  the  known  inlet  temperature.  Next,  the  same  procedure  as  in  the 
Isothermal  Model  is  used  to  calculate  S)gx£t,  and  Vg  at 

the  exit,  and  local  nozzle  area.  Finally,  the  gas  temperature  at  the  segment 
exit  is  calculated  with  the  heat  balance  of  Eq.  1.50. 

The  isentropic  efficiency  is  a  measure  of  nozzle  performance  and  is 
defined  as: 
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This  efficiency  is  the  ratio  of  the  kinetic  energy  per  unit  mass  of  the  exit 
mixture  to  the  ideal  change  of  enthalpy  (into  kinetic  energy)  for  the  isentropic 
expansion  of  a  homogeneous  two-phase  mixture.  The  maximum  enthalpy  change 
through  the  nozzle,  Ah)max,  is  given  by  the  expression: 
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The  (AP/pt)  term  in  the  second  brackets  is  included  in  the  liquid 
enthalpy  contribution  because  the  process  is  not  at  constant  (or  even  nearly 
constant)  pressure.  Although  the  temperatures  of  the  gas  and  of  the  liquid 
could  differ  at  the  nozzle  inlet,  there  is  only  one  exit  temperature  for  the 
homogeneous  mixture,  Tgx^t.  This  exit  temperature  is  calculated  by  using  the 
expression  for  an  isentropic  expansion: 
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where  R  is  the  gas  constant  and  C  is  the  mass-weighted  specific  heat: 
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(1.57) 


and  where  Tnix^in  *s  t*ie  mass-weighted  fixture  temperature  at  the 
nozzle  inlet: 
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Equation  1.56  is  derived  from  the  thermodynamic  relations  for  an  isentropic 
process  (Ref.  25),  where  the  integral  form  of  the  relationship  for  the 
entropy  change  of  a  perfect  gas  is: 
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The  corresponding  expression  for  the  entropy  change  of  a  liquid  (because 
there  is  no  change  in  volume)  is: 
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Since  the  overall  change  in  entropy  for  the  process  is  zero: 
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Therefore,  Eqns.  1.59,  1.60,  and  1.61  can  be  combined,  and  the  result 
is  Eq.  1.56,  after  the  anti-logarithm  is  taken. 

Calculated  results  using  the  Model  with  Droplet  Heat  Transfer 
for  selected  operating  conditions  are  presented  and  discussed  later  in 
this  section. 

1.4.3  Comparisons  Between  Models  and  Test  Data 

To  use  the  two  models  developed  in  Sections  1.4.1  and  1.4.2  for 
the  calculation  of  nozzle  designs,  experimental  data  are  required  in  three 
areas:  droplet  breakup  (discussed  in  Section  I. 4. 1.4  in  terms  of  the  Weber 
Number  breakup  criterion);  the  droplet  drag  coefficient;  and  the  droplet  heat 
transfer  coefficient  (required  only  for  the  Model  with  Droplet  Heat  Transfer). 

The  empirical  data  base  for  the  drag  coefficient  and  the  heat  transfer  coefficient 
have  been  developed  from  previous  investigations  published  in  the  literature. 

1 . 4 . 3 . 1  Drag_and  Heat_  Tr an j[f e^"_D at a_ 

Experimental  data  are  required  to  estimate  the  drag  coefficient  (CD) 
of  the  droplets  in  the  gas  flow.  The  drag  coefficient  describes  the  viscous 
interaction  between  the  gas  and  the  droplets  when  the  droplets  are  not  in 
velocity  equilibrium  with  the  gas.  For  solid  spheres,  Stonecypher  (as  cited 
in  Ref.  8)  made  the  following  least-squares  fit  of  the  drag  data  for  spherical 
particles  as  tabulated  by  Perry  et  al . ,  (Ref.  26): 

2  3 

in  C  =  3.271  -  0.8893  L„  +  0.03417  L  +  0.001443  L 

D  Re  Re  Re  (1.62) 

where  LRg  =  in  ^ReRej).  This  correlation  is  valid  in  the  range  of  relative 
Reynolds  Number,  0.1  <  ReRgl  <  2  x  10^  and  is  shown  in  Fig.  1.20.  Actually, 
this  correlation  applies  to  particles  at  low  concentration  in  the  gas  so  that 
no  particle  interactions  occur  and  to  flows  with  a  large  cross-section  relative 
to  the  droplet  diameter  so  that  boundary  effects  are  minimal.  If  the  loading 
ratio  is  less  than  unity,  then  the  average  particle  spacing  is  ten  times  or 
more  the  particle  size,  and  the  assumption  of  no  particle  interactions  would 
seem  reasonable  for  these  conditions. 

Ingebo  (Ref.  27)  measured  drag  coefficients  of  a  particle  cloud  injected 
into  a  wind  tunnel  at  particle  mass  fractions  of  about  one  percent,  in  the 
range  of  relative  Reynolds  Numbers  of  approximately  6  to  400.  The  data,  which 
are  correlated  by  the  expression: 
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falls  below  those  of  Stonecypher * s  least-squares  fit  as  seen  in  Fig.  1.20. 
Ingebo  attributed  the  lower  drag  coefficients  to  the  fact  that  the  particles 
were  being  accelerated,  thereby  implying  non-steady-state  flow. 

In  the  flow  behind  weak  shock  waves,  Rudinger  (Ref.  28)  measured  drag 
coefficients  for  flows  of  suspensions  of  spherical  glass  particles  in  air. 
Mixtures  of  particles  with  average  diameters  of  29  microns  (0.0011  in.)  and  of 
62  microns  (0.0024  in.)  were  tested,  over  a  range  of  loading  ratios  from  0.05 
to  0.36  and  a  range  of  relative  Reynolds  Numbers  from  12  to  600.  The  data  are 
higher  than  predicted  by  a  correlation  such  as  that  of  Eq.  1.62  (which  is  based 
on  single-sphere  data),  and  the  slope  of  the  correlation  (valid  for  50  < 

ReRej  <  300)  is  much  steeper: 
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(1.64) 


Rudinger  developed  a  model  to  account  for  these  higher  drag  coefficients 
which  were  measured  at  loading  ratios  higher  than  those  of  Ingebo.  This 
preliminary  model  allowed  for  both  microscopic  longitudinal  and  lateral  per¬ 
turbations  of  particle  motion  resulting  from  such  effects  as:  particle  inter¬ 
actions  with  the  wakes  of  other  particles;  lateral  effects  caused  by  particle 
rotation;  or  electrostatic  forces.  The  drag  coefficient  is  evaluated  for  the 
local  vectorial  particle  velocity  relative  to  the  gas  by  a  correlation  similar 
for  that  in  Eq.  1.62,  and  Rudinger  attributed  the  higher  effective  drag  co¬ 
efficients  to  these  perturbations  in  particle  motion.  In  other  words,  the 
effective  drag  coefficient  for  the  average  motion  differs  from  that  applicable 
to  the  local  microscopic  motion. 

A  rigorous  analytical  study  by  Gluckman  et  al.  (Ref.  29)  examined  the 
drag  coefficient  of  spheres  in  an  unconnected  chain  (with  the  flow  direction 
being  the  same  direction  as  the  axis  of  the  chain)  in  unbounded  Stokes  flow. 
Stokes  flow  generally  exists  at  relative  Reynolds  Numbers  below  unity.  Ana¬ 
lytically,  their  work  indicated  a  reduction  in  the  drag  coefficient  on  spheres 
in  the  middle  of  the  chain.  At  a  sphere  spacing  of  four  diameters,  the  drag  is 
about  38  percent  of  the  drag  predicted  by  the  classical  Stokes  drag  (Cp  * 
24/ReRe^);  at  a  spacing  of  two  diameters,  the  drag  is  about  24  percent;  while 
at  a  spacing  of  one  diameter  (spheres  touching),  the  drag  is  only  about  14 
percent.  Hence,  significant  reductions  in  the  drag  are  predicted,  because 
the  wake  from  one  sphere  influences  the  flow  around  the  succeeding  sphere. 
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Shown  in  Fig.  1.20  are  the  data  correlations  and  prediction  for  the 
drag  coefficient  that  were  just  discussed.  Based  on  these  data,  a  wide  range 
in  the  value  of  the  drag  coefficient  would  be  expected.  For  the  anticipated 
range  of  operating  conditions  (100  <  ReRg^  <  10,000)  of  the  two-phase 
nozzle,  the  data  are  in  the  vicinity  of  Stonecypher's  curvefit,  and  accordingly 
this  is  the  comprehensive  expression  that  is  used  in  the  computerized  models 
to  calculate  the  droplet  drag  coefficient.  However,  the  sensitivity  of  the 
nozzle  performance  predictions  to  this  assumption  for  the  drag  coefficient 
will  be  examined  in  Phase  II  of  this  study. 

Similarly,  experimental  data  are  required  to  estimate  the  heat  transfer 
coefficient  (h)  between  the  liquid  droplets  and  the  gas.  McAdams  (Ref.  19) 
presents  data  for  heat  transfer  between  spheres  and  a  gas  (air).  In  the 
range  of  relative  Reynolds  number  (  s  p  VgiipD/Pg)  25  to  10,000,  the 

dimensionless  correlation  recommended  by  McAdams  is: 
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The  relative  Reynolds  Number  for  some  of  the  nozzle  flow  conditions  of 
interest  for  marine  propulsion  will  be  in  the  range  of  100  to  2000.  For 
conditions  where  the  Reynolds  Number  is  much  lower  (e.g.,  in  the  range  of  1.0 
to  25.),  Kreith  (Ref.  30)  recommends  the  following  correlation: 
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Equations  1.65  and  1.66  are  the  two  correlations  that  are  used  to  calculate  the 
heat  transfer  coefficient  in  the  Model  with  Droplet  Heat  Transfer. 


1 . 4 . 3 . 2  £hennoj>h£S  ic  al_P£0£e  rt  i_e  s_ 

Thermo physical  properties  of  both  the  liquid  and  the  gas  are  required 
to  obtain  numerical  results  from  the  models.  Some  testing  (Ref.  8)  has  been 
done  with  mixtures  of  nitrogen  and  water  at  room  temperature,  and  the  thermo¬ 
physical  properties  of  these  fluids  at  these  conditions  are  known.  As 
discussed  in  Section  1.2.2,  a  liquid/gas  combination  of  interest  for  the 
marine  propulsion  application  is  Dowtherm  A/steam. 

The  thermophysical  properties  of  steam  (pg,  kg,  Pg,  Cpg)  are  known 
(for  example,  see  Ref.  31)  for  the  temperature  and  pressure  range  of  interest 
for  a  marine  propulsion  application.  Certain  of  the  thermophysical  properties 
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of  Dowtherm  A  (DTA)  that  are  required  in  the  model  (Pt,  C)  are  also 
known  (Ref.  32)  for  the  temperature  range  of  interest,  but  the  surface 
tension  (°)  has  apparently  been  measured  only  at  relatively  low  temperatures, 
that  is,  in  the  range  of  20  C  (68  F)  to  60  C  (140  F) .  (The  surface  tension 
of  the  liquid  is  required  to  evaluate  the  Weber  Number  droplet  breakup 
criterion.)  The  operating  temperatures  for  the  nozzle  design,  however,  are 
expected  to  be  as  high  as  316  C  (600  F) .  Hence,  the  surface  tension  data  had 
to  be  extrapolated.  Generally,  as  the  temperature  of  a  liquid  approaches  its 
critical  temperature  (which  is  498  C  (927  F)  for  DTA),  the  surface  tension 
approaches  zero,  because  the  interface  between  the  liquid  and  the  gas  disappears 
at  the  critical  point  of  the  substance.  Fig.  1.21  shows  surface  tension  data 
for  two  non-viscous  fluids,  as  presented  in  Ref.  33.  If  a  straight  line  were 
placed  through  the  data  for  each  fluid,  then  the  temperature  at  which  the 
surface  tension  becomes  zero  would  be  approximately  19  C  (35  F)  lower  than 
the  critical  temperature  of  each  substance.  As  seen  in  Fig.  1.21,  the 
low-temperature  data  for  the  surface  tension  of  DTA  was  extrapolated  by 
placing  a  line  through  the  three  low-temperature  points  and  through  the  point 
of  zero  surface  tension  at  19C  (35  F)  less  than  its  critical  temperature. 

The  sensitivity  of  the  nozzle  performance  predictions  to  this  assumption  for 
the  surface  tension  (and  indeed  to  the  assumption  of  a  critical  Weber  Number 
of  six  for  droplet  breakup)  will  be  examined  during  Phase  II  of  this  study. 

1.4. 3. 3  Domains 

As  noted  previously,  the  main  difference  between  the"  Isothermal  Model 
and  the  Model  with  Droplet  Heat  Transfer  is  that  the  latter  includes  the 
effects  of  heat  transfer,  since  both  models  allow  for  velocity  slip,  droplet 
drag,  and  droplet  breakup.  Despite  the  fact  that  conceptually,  the  Isothermal 
Model  is  relatively  simpler,  both  models  must  still  be  solved  numerically 
with  the  computer,  because  the  calculations  are  too  tedious  to  be  done  by 
hand  for  the  many  small  segments  that  are  required  for  numerical  stability. 
Therefore,  in  terms  of  obtaining  numerical  results,  the  Isothermal  Model 
offers  no  advantage  over  the  Model  with  Droplet  Heat  Transfer  because  computer¬ 
ized  solutions  are  required  in  both  instances. 

For  quick  calculations  which  can  be  done  manually,  however,  the  isothermal 
model  proposed  by  Rudinger  (Ref.  4)  and  discussed  in  Section  1.3.2  is  very 
useful  and  convenient.  In  such  a  case,  the  simplification  that  results  from 
the  assumption  of  isothermal  conditions  is  significant  enough  to  allow  hand 
calculations,  whereas  the  computer  would  be  necessary  for  solutions  without 
this  assumption.  Generally,  the  isothermal  model  provides  reasonable  results 
for  operating  conditions  of  high  loading  ratio  and  small  droplet  size.  As 
deduced  from  the  Weber  Number  criterion,  small  droplets  generally  occur  with 


44 


R80-9 54624-2 


fluids  Chat  have  low  surface  tension  and/or  at  operating  conditions  where  the 
velocity  slip  is  large,  thereby  creating  large  aerodynamic  forces  for  droplet 
breakup.  The  implication  of  a  large  loading  ratio  (r)  is  that  the  liquid 
acts  as  a  heat  source  to  the  expanding  gas,  thereby  reducing  its  temperature 
drop. 


In  contrast  to  isothermal  models  (either  Rudinger's  model  or  the 
computerized  model  developed  in  Section  1.4.1),  the  Model  with  Droplet  Heat 
Transfer  applies  to  the  full  range  of  nozzle  operating  conditions,  as  long  as 
droplet  flow  (in  contrast  to  slug  flow,  for  example)  exists  and  the  droplets  are 
uniformly  distributed  across  the  nozzle  cross-section.  Fig.  1.22  shows 
calculated  results  from  the  two  computerized  models,  for  a  mixture  of  steam 
and  Dowtherm  A  (DTA).  Here,  the  gas  and  liquid  velocities,  the  nozzle  area, 
and  the  droplet  size  are  calculated  as  a  function  of  axial  distance  through 
the  nozzle  for  the  specified  inlet  conditions  and  pressure  profile.  For 
these  operating  conditions,  the  velocity  slip  between  the  gas  and  the  liquid 
(V  -  V^)  increases  as  the  flow  proceeds  through  the  nozzle.  For  droplet 
breakup,  a  Weber  Number  equal  to  six  was  used  to  calculate  the  local  maximum 
droplet  size  throughout  the  nozzle.  The  Weber  breakup  criterion  is  not 
exact,  and  the  diameter  history  of  a  given  droplet  varying  slowly  over  time 
(Fig.  1.22)  is  an  idealization.  In  reality,  the  droplet  probably  breaks 
apart  into  smaller  droplets  at  irregular  intervals,  with  step  changes  in 
droplet  size  occurring.  At  this  relatively  high  value  of  the  loading  ratio 
(r  =  40),  the  results  from  the  two  models  for  velocity  and  throat  area 
virtually  overlap.  The  temperature  drop  of  the  liquid  is  about  5  C  (9  F), 
and  that  of  the  gas  is  about  15  C  (27  F) .  For  steam  alone  expanded  through 
this  pressure,  the  isentropic  temperature  drop  would  be  about  233  C  (420  F), 
so  that  the  liquid  DTA  does  act  as  a  heat  source  for  the  steam. 

Table  1.4  compares  results  calculated  from  the  Isothermal  Model  and  from  the 
Model  with  Droplet  Heat  Transfer,  at  the  same  inlet  conditions  (temperature 
and  pressure)  and  the  same  exit  pressure,  for  four  values  of  the  loading 
ratio.  At  a  given  loading  ratio,  the  exit  velocities  calculated  with  the 
Heat  Transfer  Model  are  lower  than  those  calculated  with  the  Isothermal 
Model.  Understandably,  both  models  are  idealized,  but  the  model  with  Droplet 
Heat  Transfer  is  more  realistic  than  the  Isothermal  Model  because  the  effect 
of  finite  heat  transfer  is  included  in  the  model  with  Droplet  Heat  Transfer. 

The  inclusion  of  this  heat  transfer  effect  means  that  the  calculated  velocities 
and  efficiency  are  more  realistic  and  lower  than  those  calculated  with  the 
Isothermal  Model,  which  actually  requires  heat  additions  to  the  nozzle  to 
maintain  isothermal  conditions.  Because  both  the  gas  and  the  liquid  are 
calculated  as  being  lower  in  temperature  in  the  Model  with  Droplet  Heat 
Transfer,  the  gas  is  denser.  For  a  fixed  pressure  drop  (which  translates 
into  a  fixed  amount  of  momentum),  the  gas  velocity  would  logically  be  ex¬ 
pected  to  be  lower  when  the  gas  density  is  larger,  and  the  liquid  velocity 
would  also  be  lower  because  it  is  related  to  the  gas  velocity  by  the  droplet 
drag  relationship. 
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With  the  Model  with  Droplet  Heat  Transfer,  it  can  be  seen  in  Table  1.4 
that  when  the  loading  ratio  is  increased,  both  the  gas  and  liquid  exit 
velocities  decrease  while  the  exit  temperatures  for  both  increase.  As  the 
loading  ratio  increases,  the  liquid  phase  increasingly  carries  more  of  the 
total  momentum,  but  with  the  much  higher  density  of  the  liquid,  both  velocities 
can  be  lower  while  the  mixture  still  has  the  same  total  momentum.  Similarly, 
as  the  loading  ratio  is  increased,  the  heat  capacity  of  the  liquid  phase  is 
larger,  and  the  gas  temperature  decrease  during  the  expansion  is  counteracted 
by  heat  transferred  from  the  larger  heat  source  of  the  liquid  phase. 

The  isentropic  nozzle  efficiencies  (n)  as  calculated  with  the  Model  with 
Droplet  Heat  Transfer  are  lower  than  those  with  the  Isothermal  Model.  The 
isentropic  nozzle  efficiency  used  here  is  the  ratio  of  the  kinetic  energy  per 
unit  mass  of  the  exit  mixture  to  the  ideal  change  of  enthalpy  (into  kinetic 
energy)  for  the  isentropic  expansion  of  a  homogeneous  mixture.  Because  the 
efficiency  is  a  measure  of  how  much  pressure  energy  is  converted  into  kinetic 
energy  (velocity),  the  lower  efficiency  is  explained  by  the  realistically 
lower  velocities  calculated  using  the  model  with  Droplet  Heat  Transfer. 

1.4. 3.4  Comparison  With  Other  Mode_l_s_and_Te_st^  I)a£a_ 

The  results  from  models  developed  during  this  program  can  be  compared 
with  results  of  models  previously  developed  by  other  investigators  and  with 
pertinent  experimental  data.  Table  1.5  summarizes  selected  results,  at  one 
set  of  operating  conditions,  from  the  computerized  Isothermal  Model  described 
in  Section  1.4.1,  from  Rudinger's  analytical  isothermal  model  (Ref.  4),  and 
from  the  isentropic  model  discussed  in  Section  1.3.1.  The  velocities  calcu¬ 
lated  with  the  isentropic  model  are  lower  in  all  cases  because  the  mixture 
temperature  is  lower,  hence  the  gas  is  denser,  and  the  velocity  can  be  lower 
and  still  carry  the  all-liquid  momentum.  For  the  all-gas  case  (m^/m^  *  0) 
and  the  nearly  all-liquid  case  (mj/m^  ■  1000),  the  results  for  exit 
velocity  from  the  computerized  model  and  from  Rudinger's  model  are  in  close 
agreement.  Figure  1.23  shows  the  all-gas  results  for  velocity  and  area  as  a 
function  of  distance  through  the  nozzle  for  the  two  isothermal  models;  the 
calculated  results  agree  within  one  percent.  Essentially,  this  comparison  of 
the  all-gas  results  from  the  two  models  represents  a  check  on  the  formulation 
of  the  momentum  equation  in  the  computerized  Isothermal  Model.  For  the  cases 
with  mj/m^  *  40,  the  difference  between  the  exit  velocity  results  can 
be  attributed  to  the  varying  V^/V  ratio  in  the  computerized  model,  when 
compared  to  the  condition  that  V^/V  *  constant  in  Rudinger's  model.  This 
assumption  of  V^/Vg  “  constant  replaces  the  droplet  drag  equation  that  couples 
the  liquid  and  the  gas  momentums  in  the  computerized  model.  Further,  Rudinger's 
model  is  for  constant-diameter  (solid)  particles. 
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Table  1.6  compares  velocity  and  temperature  results  at  the  exit  of 
a  50-inch-long  nozzle  (calculated  by  Elliot  and  Weinberg  (Ref.  8)  at  their 
selected  conditions)  with  the  results  from  three  other  models  for  the  two-phase 
flow  of  a  mixture  of  water  and  nitrogen  at  inlet  conditions  of  150  psia  and 
16  C  (60  F) .  The  results  from  Rudinger's  isothermal  model  are  shown  at 
V^/Vg  =  0.7,  the  value  at  the  nozzle  exit  which  is  also  the  approximate 
value  throughout  the  nozzle  obtained  from  the  results  of  Elliott  and  Weinberg. 
The  exit  velocity  (V)  calculated  with  either  of  the  two  isothermal  models 
is  higher  than  the  exit  velocity  calculated  with  either  of  the  other  two 
models,  and  hence,  the  isentropic  nozzle  efficiencies  calculated  with  the  tw^ 
isothermal  models  are  higher  also.  As  discussed  earlier,  both  the  Isothermal 
Model  and  the  Model  with  Droplet  Heat  Transfer  are  best  solved  by  computer, 
but  Rudinger's  isothermal  model  is  suitable  for  rapid  manual  calculations, 
although  with  the  latter,  the  calculated  nozzle  efficiencies  will  be  opti¬ 
mistically  high. 

Results  from  the  Model  with  Droplet  Heat  Transfer  are  shown  in  Table  1.6 
for  three  different  pressure  profiles  selected.  The  profile  with  dP/dX  * 
constant  is  linear  with  distance  between  the  inlet  and  the  exit,  while  the 
2-line  and  3-line  curvefits  are  shown  in  Figure  1.24  along  with  the  optimum 
pressure  profile  for  maximum  mean  exit  velocity  as  calculated  by  Weinberg  and 
Elliott.  The  purpose  behind  using  these  selected  pressure  profiles  is  to 
compare  the  results  calculated  from  the  Model  with  Droplet  Heat  Transfer  with 
the  results  presented  by  Weinberg  and  Elliot,  which  correspond  to  the  optimum 
pressure  profile.  The  nozzle  cross-sectional  area  profiles  which  correspond 
to  each  of  the  three  pressure  profiles  shown  in  Fig.  1.24  and  to  the  linear 
pressure  profile  are  all  slightly  different.  The  liquid  exit  velocity  (and 
the  weighted  mean  velocity),  the  exit  temperature  of  the  liquid  and  of  the 
gas,  and  the  isentropic  nozzle  efficiency,  all  of  which  are  calculated  with 
either  of  the  three  pressure  profiles  in  the  computerized  Model  with  Droplet 
Heat  Transfer  or  with  Elliott  and  Weinberg's  model  and  their  optimum  pressure 
profile,  agree  within  one  percent.  Furthermore,  the  liquid  carries  most  of 
the  momentum  and  enthalpy  at  the  loading  ratios  expected  for  the  two-phase 
nozzle  in  a  possible  marine  propulsion  application.  The  gas  velocities, 
which  were  calculated  in  the  same  manner  with  the  two  models,  differ  by  18 
percent.  Apparently,  the  gas  velocity  results  are  much  more  sensitive  to  the 
pressure  profile  selected  than  are  the  liquid  velocity  results,  as  the 
results  for  the  gas  velocity  calculated  with  the  computerized  Isothermal 
Model  (which  is  just  the  momentum  equation  and  no  heat  transfer  equation)  and 
seen  in  Table  1.5  for  the  same  linear  pressure  profile,  are  close  to  the 
results  from  the  analytical  isothermal  model  of  Rudinger. 

Both  the  Model  with  Droplet  Heat  Transfer  and  the  model  developed  by 
Elliott  and  Weinberg  use  the  basic  principles  of  continuity,  momentum,  droplet 
drag,  droplet  breakup,  energy,  and  heat  transfer.  However,  the  formulation 
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of  the  Model  with  Droplet  Heat  Transfer  emphasizes  its  suitability  to 
calculate  parametric  results,  whereas  the  formulation  of  the  Elliott  and 
Weinberg  model  emphasizes  nozzle  design  details  but  is  too  complicated  for 
parametric  analyses  when  many  design  cases  need  to  be  considered  in  a  short 
period  of  time.  The  Elliott  and  Weinberg  design  procedure  consists  of 
several  steps.  First  the  optimum  nozzle  contour  is  determined  with  their 
computerized  model  at  a  condition  of  constant  droplet  diameter  (because  the 
pressure  profile  may  be  so  steep  at  the  nozzle  inlet  as  to  prevent  droplet 
breakup).  Computer  runs  are  made  for  a  series  of  values  of  the  optimizing 
Lagrangian  multiplier,  X  ,  to  cover  the  range  of  nozzle  lengths  of  possible 
interest.  Second,  a  preliminary  guess  of  the  optimum  nozzle  length  is  made 
to  permit  modifying  the  optimum  cross-sectional  profile  into  a  practical 
nozzle  contour.  Nozzles  that  are  too  short  may  have  large  velocity  slip 
while  nozzles  that  are  too  long  may  have  large  friction  losses;  the  result  in 
either  extreme  case  is  a  low  mean  exit  velocity  (and,  hence  low  nozzle 
efficiency).  This  procedure,  while  eminently  suitable  for  nozzle  design,  is 
too  detailed  for  evaluating  numerous  possible  nozzle  operating  conditions. 

In  contrast,  the  Model  with  Droplet  Heat  Transfer  is  formulated  with 
regard  to  parametric  analysis,  but  local  velocity  and  temperature  conditions 
are  also  calculated  and  available  throughout  the  nozzle.  This  model  neglects 
friction,  but  this  assumption  is  reasonable  since  the  nozzles  for  a  two-phase 
turbine  application  must  be  short  in  length  in  order  to  accommodate  them 
into  the  available  space.  Further,  frictional  effects  did  not  appear  to  be 
large  even  in  results  for  the  50-inch-long  nozzle  of  Table  1.6. 

Experimental  data,  as  discussed  in  Section  1.1.1,  have  been  taken  for 
a  number  of  two-phase  nozzles.  However,  only  some  of  these  data  were  taken 
with  a  two-component  mixture,  the  type  of  system  to  which  the  Model  with 
Droplet  Heat  Transfer  applies.  A  complication  arises  when  comparing  experi¬ 
mental  data  for  the  performance  (both  design-point  and  off-design)  of  a  given 
nozzle  with  the  calculated  results  from  a  design  model.  Essentially,  a 
design  model  is  used  to  plan  families  of  nozzles  for  one  set  of  operating 
conditions,  whereas  a  performance  model  is  used  in  the  reverse  manner,  that 
is,  to  calculate  the  performance  of  one  nozzle  at  many  sets  of  operating 
conditons.  Experimental  performance  data  (Ref.  7)  for  a  six-inch  nozzle 
tested  with  a  mixture  of  water  and  air  are  shown  in  Fig.  1.25.  Also  shown 
are  performance  lines  predicted  for  the  actual  nozzle  contour  by  Elliott  and 
Weinberg  with  their  model  (Ref.  B).  The  Elliott  and  Weinberg  prediction 
labeled  "Real"  is  based  on  a  model  that  includes  velocity  slip,  droplet 
breakup,  and  heat  transfer,  but  no  friction;  while  their  prediction  labeled 
"With  Friction"  includes  friction  in  terms  of  wall  boundary  layer  losses. 
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Both  of  the  computerized  models  presented  in  Sections  1.4.1  and  1.4.2 
were  also  used  to  make  nozzle  designs  for  a  range  of  loading  ratios  at  the 
inlet  temperatures,  velocities,  and  pressure  of  this  six-inch  nozzle,  by 
assuming  a  linear  variation  of  pressure  drop  with  distance  through  the 
nozzle.  (As  discussed,  the  nozzle  contour  produced  from  these  design  models 
would  be  different  at  each  value  of  the  loading  ratio.)  The  results  from 
these  design  models  are  also  shown  in  Fig.  1.25.  The  design-point  performances 
of  a  series  of  nozzles  compare  well  with  the  experimental  data  and  show  the 
same  trend  as  the  off-design  performance  of  one  nozzle.  The  predictions  from 
the  Model  with  Droplet  Heat  Transfer  are  realistically  lower  than  those  from 
the  Isothermal  Model,  and  the  difference  is  larger  at  lower  values  of  the 
loading  ratio  where  the  temperature  of  the  gas  phase  after  expansion  is  lower 
because  the  gas  represents  a  relatively  larger  heat  sink  at  lower  values  of 
the  loading  ratio. 

1.4.4  Preliminary  Parametric  Calculations 


As  part  of  an  effort  to  check  the  computer  models,  some  preliminary 
parametric  calculations  were  made  to  show  their  versatility.  Figure  1.26 
shows  the  effect  of  inlet  temperature  on  temperatures,  velocities,  and 
cross-sectional  area  through  the  nozzle.  In  one  case,  the  inlet  temperatures 
of  the  gas  and  the  liquid  are  the  same  at  316  C  (1060  R) ,  while  in  the  other 
case  the  inlet  gas  temperature  is  56  C  (100  R)  cooler  than  the  inlet  liquid 
temperature  of  1060  R.  The  liquid  inlet  temperatures  were  selected  to  be  the 
same  because  the  liquid  contains  the  bulk  of  the  mixture  enthalpy  at  a 
loading  ratio  of  40.  Near  the  inlet  for  the  case  with  the  100  R  differential, 
the  steam  is  rapidly  heated  by  the  higher-temperature  liquid  DTA,  until  the 
steam  expands  sufficiently  for  both  temperatures  to  begin  falling  shortly 
beyond  the  nozzle  exit.  The  results  for  the  liquid  velocity  in  both  cases 
overlap,  as  also  do  the  gas  velocities.  For  the  case  with  the  cooler  gas, 
slightly  less  nozzle  area  is  required,  because  the  cooler  gas  is  denser  than 
in  the  case  with  the  higher  inlet  temperature. 

The  effect  of  nozzle  length  on  performance  is  shown  in  Table  1.7.  As  the 
length  is  increased,  the  residence  time  of  the  two-phase  mixture  in  the  nozzle 
increases,  and  more  enthalpy  and  momentum  can  be  transferred  between  the  gas 
phase  and  the  liquid  phase.  The  weighted  exit  velocity  is  higher,  and  as  a 
result,  the  nozzle  efficiency  increases  with  increased  nozzle  length.  At 
this  value  of  the  loading  ratio,  the  liquid  temperature  decreases  only  0.3  C 
(0.6  R)  as  the  nozzle  length  is  increased  from  3  in.  to  12  in;  for  these 
operating  conditions,  the  nozzle  efficiency  increases  by  three  percentage 
points  with  the  added  length. 
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The  effect  of  loading  ratio  (m^/nig)  on  nozzle  performance 
is  shown  in  Table  1.4  for  one  set  of  operating  conditions.  As  the  loading 
ratio  increases,  the  two-phase  flow  in  the  nozzle  tends  to  become  increasingly 
isothermal,  because  the  thermal  capacity  of  the  liquid  phase  dominates. 

The  preliminary  parametric  calculations  presented  here  indicate  that 
the  computerized  Model  with  Droplet  Heat  Transfer  is  appropriately  formulated 
for  parametric  analyses  of  two-phase  nozzles. 
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SUMMARY  OF  ANALYTICAL  TWO-PHASE-NOZZLE  MODELS 
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TABLE  1.2 

POTENTIAL  TWO-PHASE  NOZZLE  APPLICATIONS 

Low-RPM  Two-Phase  Turbines 
Marine  Propulsion 
Automobile 

Two-Phase  Turbines  for  Geothermal  Power 
Water  Injection  for  Jet  Engine  Thrust  Augmentation 
Two-Phase  Nozzles  for  Liquid-Metal  MHD  Power 
Two-Phase  Turbines  for  Diesel  Bottoming 
Two-Phase  Water  Piston  Turbines 
Rocket  Thrustors 

Gas-Driven  Jet  Pumps  for  Rocket  Engines 
Two-Phase  Water  Jet  Thrustors 
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TABLE  1.3 

PARAMETERS  OF  PREVIOUS  TWO-PHASE-NOZZLE  INVESTIGATIONS 


Investigator 

Gas/Liquid 
Combinat ion 

Inlet 

Temp/Press 

Liquid 

Mass  Fraction 

Liquid 

Vol.  Fraction 

Droplet 

Diameter 

Alger 

Geothermal 

Steam/Water 

435  F/360 
psia 

0.81  inlet 

0.06  inlet 

1  to  10  urn 

Comfort  & 

Beadle 

Geothermal 

Steam/Water 

434  F/360 
psia 

0.81  inlet 

0.65  exit 

0.06  inlet 
0.01  exit 

6.5  urn 

Elliot  & 

Hayes 

Power  Turbine 
Steam/oil 

620  F/1450 
ps  ia 

0.95 

+  25  um 

Elliot  & 

We inberg 

MHD  Power 
Cesium/Lith ium 

1800  F/137 
psia 

0.93 

0.28 

1270  um 

Nitrogen/Water 

60  F/150 
psia 

0.976 

0.33 

780  Um 

Freon  1301/ 
Water 

50  F/150 
psia 

0.75  to 

0.93 

Elliot 

Air/Water 

50  F/514 
psia 

0.75  to 

0.97 

Starkman 

et  al 

Steam/Water 

580  F/1000 
ps  ia 

0.8  to 

1.0 

250  um 

Schrock 

et  al 

Steam/Water 

550  F/1305 
ps  ia 

0.9  to 

1.0 

Ahmad  & 
Hayes 


Bottoming  Cycle  660  f/2365 
Steam/Dowtherm  psia 
A 


0.8 
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Table  1.4 

DROPLET  HEAT  TRANSFER  MODEL  VERSUS  ISOTHERMAL  MODEL 


Six-Inch  Nozzle 


Steam  ^)TA  Inlet:  1060  R  and  1000  psia 

Exit:  100  psia 


At  Exit 


Model 

m2./mg 

vs 

v* 

V 

Tg 

T* 

n 

ft/sec 

- 

°R 

— 

Heat  Transfer 

5 

1665 

1383 

1430 

957 

1009 

0.89 

Isothermal 

5 

1698 

1417 

1464 

— 

1060 

0.93 

Heat  Transfer 

10 

1312 

1065 

1087 

993 

1030 

0.88 

Isothermal 

10 

1326 

1080 

1102 

— 

1060 

0.90 

Heat  Transfer 

20 

1036 

822 

832 

1017 

1043 

0.88 

Isothermal 

20 

1041 

829 

839 

1060 

0.89 

Heat  Transfer 

40 

838 

654 

658 

1033 

1051 

0.90 

Isothermal 

40 

840 

656 

661 

1060 

0.91 
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Table  1.5 


COMPARISON  OF  ISOTHERMAL  MODELS  (Tj;,=Tg=Tinlet=constant) 


6-inch  Nozzle 
Steam/DTA 


Model 


/  Wg 


Isothermal 

Rudinger 

Isentroplc 

Isothermal  (Varying  V^/V  ) 
Rudinger  (Vj/Vg=const=0. B3) 
Rudinger  (V£/Vg=const=0. 78) 
Isentropic 

Isothermal 

Rudinger  (Vj=Vg=const=0. 82) 
Isentropic 


Inlet:  1060  R,  1000  psia 
Exit:  100  psia 


■At  Exit- 


vg 

V 

ft/sec 

0 

0 

3674 

- 

3674 

0 

0 

3674 

- 

3674 

0 

0 

3245 

- 

3245 

40 

0.78 

840 

656 

661 

40 

0.83 

792 

657 

660 

40 

0.78 

824 

643 

657 

40 

1.00 

571 

571 

571 

1000 

0.82 

508 

417 

417 

1000 

0.82 

510 

418 

418 

1000 

1.00 

116 

116 

116 

>A11-Gas 


^Nearly 
f All-Liquid 
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Table  1.6 

COMPARISON  OF  FOUR  MODELS 


50-inch  Nozzle  with  Nitrogen/Water 

Inlet:  520  R  and  150  psia 
Exit:  14.1  psia 

«\/mg  =  40 


Exit  Values 


Model 

Vg 

V£ 

V 

Tg 

B 

B 

ft/sec 

0 

R 

Rudinger  Isothermal  Model 

VVg  -  0.7 

437 

306 

— 

309 

520 

0.758 

Isothermal  Model 

Varying  V£/Vg 

502 

301 

306 

520 

0.743 

Droplet  Heat  Transfer  Model 
dP/dX  *  const 

496 

300 

305 

503 

518.6 

0.736 

2-line  Curvefit  to 

Elliott  and  Weinberg  P  Profile 

504 

298 

302 

502 

518.6 

0.726 

3-line  Curvefit  to 

Elliott  and  Weinberg  P  Profile 

462 

300 

303 

505 

518.6 

0.730 

Elliott  and  Weinberg  Model 

Optimum  dP/dX 

425 

300 

303 

508 

518.3 

0.729 

•MW"**"  »*«- 


PHASE  MIXTURE,  T/T0 


TEMPERATURE,  T 
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FIG.  1.3 


T-S  CYCLE  PATHS  FOR 
TWO-PHASE  AND  SINGLE-PHASE  TURBINES 


78-1 


1-08-6 


WEBER  NUMBER 
DROPLET  BREAKUP 


We)MAX  _  6 


We 


vvv2° 

2  O 


STEAM/DOWTHERM  A 
600  “F 


\ 


0  8  =  V„/Vn 


VAPOR  PRESSURE  RELATIONSHIPS  FOR  VARIOUS  FLUOROCARBONS  AND  WATER 


FIG 


F 
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B!Sd  '3ynSS3dd  3imOS8V 


79-11-106 


TEMPERATURE 


PARTICLE  VOLUME  FRACTION, 


1-2 


RELATIONSHIP  BETWEEN  LOADING  RATIO,  DENSITY  RATIO, 
AND  PARTICLE  VOLUME  FRACTION 


VELOCITY 


NOZZLE  AREA,  A/A„;  PARTICLE  VOLUME  FRACTION 


GAS  EXIT  VELOCITY,  U„/Un) 
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THROAT  PRESSURE  RATIO.  P(/P0 


4-^ 


EFFECT  OF  INLET  PARTICLE  VOLUME  ON  THROAT  PRESSURE  RATIO 


VELOCITY,  tt/sec 


R80-954€ 


VELOCITY,  ft/sec 


I. 


COMPARISON  OF  HOMOGENEOUS  WITH  PARTIALLY  SEPARATED  MODEL 


78-1 


1-106 


m/sec 


R80-954624-2 


FIG.  1.14 


SIMPLIFIED  MODEL  TO  ESTIMATE  HEAT  TRANSFER  RATE  EFFECTS 


V’  OUT 

g;  out 
Tg>  out 


79-  11-106-5 


STEAM  TEMPERATURE, F 
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FIG,  1.15 


EFFECT  OF  DROPLET  SIZE  AND  VELOCITY  ON  DROPLET  HEAT  TRANSFER 


LOADING  RATIO, rfy/rfig 


STEAM  TEMPERATURE,  F 
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FIG.  1.16 


EFFECT  OF  LIQUID-GAS  TEMPERATURE  DIFFERENTIAL  ON  DROPLET  HEAT  TRANSFER 


LOADING  RATIO,  m^/rrig 


79—  1 1  —  106—3 


DIMENSIONLESS  RESIDENCE  TIME,tres//rv 


2 


FIG.  1.19 


FORCE  BALANCE  WITH  DRAG  ON  LIQUID  DROPLET 


FLOW  DIRECTION 


FDRAG—  ^ PRESS  "  FACCEL 


80-5-123-® 


COMPARISON  OF  DRAG  COEFFICIENTS 


RELATIVE  REYNOLDS  NUMBER 


Dimensionless 
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